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PREFACE 


We  became  interested  in  the  anemoclinometer  as  a 
possible  three-dimensional  pressure-sphere  anemometer  in 
1961  when  Professor  H.  Lettau  called  it  to  our  attention. 
Several  features  were  of  interest;  the  small  size  which 
would  enable  measurements  near  the  ground,  the  internal 
angular  precision  of  construction,  and  the  stability  and 
ruggedness  of  the  probe  were  all  valuable  attributes. 

Most  importantly,  however,  the  vertical  velocity,  which 
is  the  most  critical  measurement,  is  obtained  from  a 
pressure  proportional  to  the  product  of  the  vertical  and 
horizontal  winds.  Consequently,  the  vertical  wind  is  con¬ 
tained  in  a  term  of  large  magnitude,  which  can  be  measured 
with  precision,  and  is  then  obtained  by  division  rather 
than  by  differencing,  which  also  lends  precision  to  the 
measurement. 

The  early  work,  which  validated  the  potential  of  the 
anemoclinometer  as  a  three-dimensional  anemometer,  was 
supported  under  grant  DA-SIG-36-039-62-G25  by  the  Atmos¬ 
pheric  Sciences  Laboratory,  U.  S.  Army  Electronics  Command, 
Fort  Huachuca,  Arizona  (formerly  the  Department  of  Meteorol¬ 
ogy,  U .  S.  Army  Electronics  Research  and  Development  Activity, 
Fort  Huachuca,  Arizona) .  This  work,  which  included  tests 
on  frequency  response,  sensitivity  of  the  anemometer  to 
angular  rotation,  and  a  limited  comparison  to  wind  profile 
measurements  of  shear  stress,  was  reported  — ^  for  the  above 


— ^  Thurtell,  G.  W.  and  C.  B.  Tanner.  1965.  Momentum  Transport 
Measurement  in  the  Atmospheric  Surface  Layer  with  the 
Anemoclinometer.  University  of  Wisconsin,  Department  of 
Soil  Science,  Madison,  Wisconsin,  Final  Report  1962-1965. 
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grant.  The  results  also  are  presented  by  Thurtell  -  . 

The  use  of  the  pressure  probe  hinged  upon  a  pressure 
measurement  with  severe  requirements  of  sensitivity,  zero 
stability,  sensitivity  stability  ard  frequency  response. 

The  only  pressure  transducer  available  at  that  time  which 
appeared  to  meet  our  requirements  was  one  made  by  Datametrics, 
Inc.,  as  described  in  this  report  and  the  earlier  one. 
Datametrics  was  most  helpful  in  modifying  the  sensor  and 
the  electronics  to  meet  our  requirements,  including  smaller 
transducer  volume,  distant  separation  of  the  transducer 
from  the  electronics,  and  read  out  of  electrical  zero  and 
full  scale. 

Earlier  experience  with  data-logging  via  a  magnetic 
12/ 

tape  system  — convinced  us  that  the  only  feasible  route 
was  on-line  computation.  This  was  done  in  the  experiments 
discussed  in  this  report,  and  proved  to  be  as  valuable  as 
we  anticipated. 

Most  of  the  results  in  this  report  were  obtained  as 
part  of  the  1967  Cooperative  Field  Experiment  conducted 
at  the  University  of  California  at  Davis  and  sponsored  by 
the  Atmospheric  Sciences  Laboratory,  U.  S.  Army  Electronics 
Command,  Fort  Kuachuca,  Arizona.  The  remaining  data  were 
gathered  at  the  University  of  Wisconsin  Hancock  Experiment 
Farm. 

We  wish  to  express  our  appreciation  to  Mr.  T.  A.  Black# 
graduate  student  who  worked  with  us  on  the  Davis  California 
experiment  and  helped  also  at  Hancock  in  providing  all  of 


--  Thurtell,  G«  VJ.  1965.  Momentum  transport  measurements  in 
the  atmospheric  surface  layer  with  the  anemoclinometer . 
Ph.D.  Thesis,  Univ.  Wis.  (Pub.  No.  65-11,179).  47  p.  Univ. 

Microfilms,  Ann  Arbor,  Mich. (Diss.  Absfcr.  1:  4017.). 
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the  energy  balance  data.  We  owe  thanks  also  to 
Dr.  C.  R.  Stearns  (Department  of  Meteorology,  University 
of  Wisconsin) ,  Dr.  W.  0.  Pruitt  (Department  of  Water 
Science  and  Engineering,  University  of  California-Davis) , 
and  Dr.  J.  A.  Businger  (Department  of  Atmospheric  Sciei.jes, 
University  of  Washington) ,  who  participated  in  che  1967 
Cooperative  Field  Experiment  and  provided  data  used  for 
independent  comparisons  of  shear  stress  and  sensible  heat 
flux  density. 


G.  W.  Thurtell 
C.  B.  Tanner 
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THREE-DIMENSIONAL  PRESSURE-SPHERE 
ANEMOMETER  SYSTEM 


G.  W.  Thurtell,  C.  B.  Tanner,  and  M.  L.  Wesely 
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ABSTRACT 

A  rugged  and  stable  pressure-sphere  anemometer 
system  is  described  which  provides  an  accurate  measure¬ 
ment  of  wind  velocity  and  direction  within  a  meter  of 
the  ground.  The  horizontal  wind  velocity,  (u2  +  v2) 
agreed  very  closely  with  cup  anemometer  measurements, 
indicating  good  accuracy  in  the  measurement  of  the 
dominant  term,  u.  Eddy  correlation  measurements  of 
shear  stress  with  the  pressure-sphere  agreed  very  well 
with  Davis  shear-stress  meter  measurements  and  satis¬ 
factory  agreement  was  found  with  data  obtained  from 
wind  velocity  profiles  and  from  wind  measurements  using 
a  drag  coefficient.  Ratios  of  <yw/u*  during  neutral 
periods  were  found  to  be  in  excellent  agreement  with 
values  derived  by  Panofoky  and  Lettau,  providing  further 
indication  of  the  accuracy  obtainable  with  the  pressure 
sphere  system. 
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1.  Introduction 

The  basic  mechanisms  of  turbulent  transport  in  the 
layers  of  air  within  a  few  meters  of  the  earth's  surface 
are  receiving  increasing  attention  from  researchers 
from  many  disciplines.  Inadequate  diffusion  models  are 
limiting  research  progress  in  meteorology,  ecology, 
agriculture,  water  resources  and  air  pollution  since 
many  of  the  critical  problems  in  these  fields  are 
associated  with  the  exchange  of  energy,  gases  and  aer¬ 
osols  between  the  earth  and  its  atmosphere.  The  testing 
and  development  of  improved  transport  models  requires 
accurate  experimental  data  which  is  at  present  insuf¬ 
ficient. 

Field  measurements  of  turbulent  mixing  processes 
have  been  few  and  generally  inadequate  because  of  the 
stringent  requirements  for  the  instrumentation.  The 
wind  velocity  sensors  must  be  accurate,  stable  under 
field  conditions,  fast  responding,  small  for  measure¬ 
ments  near  the  ground,  rugged,  and  must  measure  both  the 
flow  direction  and  velocity  without  seriously  disturb¬ 
ing  the  flow.  Sonic  anemometers  (Kaimal,  et  al.,  1964; 
Kaimal,  et  al.,  1968),  bivanes  (Gill,  1963;  McCready 
and  Jex,  1964;  Cramer,  et  al,,  1961),  two  types  of  heat- 
transfer  anemometers;  (Miyake  and  Badgley,  1967;  Dyer, 

1960) ,  fast-response  cup  anemometers  (Frenzen,  1965)  and 
vertical,  propeller-type  anemometers  (Thorn thwa it e,  et  al., 
1961  and  Holmes,  et  al.,  1964)  have  all  been  used  in  the 
atmospheric  surface  layer  but  each  fails  to  meet  one  or 
more  of  the  essential  criteria  mentioned  above  for 
measurements  near  the  ground.  The  pressure  sphere  is 
well-suited  to  measuring  the  lateral  and  vertical  wind 
components  because  they  appear  as  products  with  the 
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large  longitudinal  velocity  in  the  basic  pressure  measure¬ 
ment.  It  is  felt  that  the  anemometer  system  to  be 
described  does  meet  these  requirements  to  greater  degree 
than  do  other  available  instruments  and  will  thus  aid 
research  on  turbulent  diffusion  processes. 

The  basic  sensor  of  our  system  is  the  anemoclincmeter 
described  by  Martinot-Lagarde,  et  al.#  (1952)  and  made 
by  the  Institut  de  Mecanique  des  Fluides-^.  The  tests  to 
be  described  were  conducted  at  the  University  of  California 
at  Davis  as  part  of  the  1967  Cooperative  Field  Experiment 
sponsored  by  the  Atmospheric  Sciences  Laboratory,  U.  S. 

Army  Electronics  Command,  Ft.  Huachuca,  Arizona.  Wind 
velocity  measurements  made  with  our  anemometer  system 
were  compared  with  cup  anemometers.  Eddy  correlation 
shear  stress  measurements  were  compared  both  with  wind 
profile  data  and  with  data  from  the  large  Davis  shear- 
stress  meter  (Brooks  and  Pruitt,  1966) .  in  addition 
measurements  of  the  standard  deviation  of  the  vertical 
component  of  wind  velocity  are  presented  for  Davis  and 
also  some  from  Hancock,  Wisconsin. 

2.  Anemometer  system 

The  anemometer  system  consists  of  a  spherical  probe 
with  pressure  ports  drilled  into  its  surface.  This 
particular  design  of  sphere,  called  Philip  I,  can  be 
replaced  satisfactorily  by  other  styles.  The  pressures 
developed  at  these  ports  are  transmitted  through  small 
tubes  and  measured  by  pressure  transducers.  The 
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Institut  de  Mecanique  des  Fluides  de  Lille,  5, 
Boulevard  Painleve,  Lille  (Nord) ,  France. 


electrical  outputs  of  the  pressure  transducers  can  be 
analyzed  to  give  the  orthogonal  components  of  the  wind 
vector. 

a .  Spherical  pressure  probe 

We  used  both  3-cir.  and  8-cm  pressure  probes  as 
described  by  Martinot-Lagarde,  et  ai.,  (1952).  The  3-cm 
prcbe,  shown  in  Pig.  1,  consists  of  a  spherical  head 
mounted  on  a  supporting  shaft.  A  drawing  of  the  head, 
showing  some  of  the  ports,  is  given  in  Fig.  2.  When  in 
use,  the  probe  is  fixed  with  the  sphere  on  the  upstream 
end  of  the  shaft.  Twelve  small  ports  are  drilled  into 
the  spherical  surface  and  a  pitot  tube  is  centered  in  a 
Venturi  which  is  bored  in  the  sphere  on  the  axis  of  the 
shaft.  Eight  of  the  twelve  ports  on  the  surface  of  the 
sphere  are  located  on  a  circle  at  an  angle  of  47.5°  to 
the  shaft  axis  and  serve  as  reference  ports  for  the 
pitot  tube  in  the  Venturi;  these  eight  holes  are  con¬ 
nected  to  a  common  pressure-averaging  cavity  in  the 
shaft.  The  other  four  holes  lie  at  right  angles  in  the 
x,z  and  x,y  planes,  and  each  hole  is  at  45°  from  the 
shaft  axis.  The  x-coordinate  is  taken  parallel  and  the 
y-  and  z- coordinates  perpendicular  to  the  probe  shaft, 
with  z  in  the  vertical  plane.  The  open  end  of  the 
pitot  tube  in  the  Venturi  is  in  the  upstream  direction. 
The  pressure  difference 


P 


1 


P 


m 


(1) 


between  the  pitot  tube  and  the  cavity  common  to  the 
eight  reference  ports,  is  proportional  to  the  dynamic 
pressure. 


Spherical  sensing  head  of  anemo- 
clinometer  showing  pressure  ports 


Front  and  cross-section  views  of 
anemoclinometer  head,  with  y-  and 
z-coordinates  shown  on  front  view 
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Pj_  =  ^apV2  (2) 

where 

2  2  2  2 

V  =  (u  +  v  +  wz)  (3) 

where  u,  v,  and  w  are  the  axial,  cross-horizontal,  and 
vertical  wind  components,  o  is  the  fluid  density,  and  "a" 
is  a  constant  of  the  probe  equal  to  1.015  according  to 
data  supplied  by  the  manufacturer.  The  pressure  differ¬ 
ence  between  the  two  vertical  ports  (x,z  plane)  and  that 
between  the  horizontal  ports  (x,y  plane)  are  predicted 
reasonably  well  by 


P2  =  bouw 


(4) 


P.^  =  bouv 


(5) 


The  factor,  b,  is  a  function  of  the  Reynold’s  number  but 
is  relatively  constant  in  the  Reynold  number  range  of 
2,000  to  200,000. 

Calibration  data  supplied  by  the  manufacturer  indi¬ 
cated  that  for  the  3-cm  spheres,  the  pressure  ratios 
P2^P1  and  P3//pl  were  linearly  related  to  the  angles  F  and 
G  respectively  by  the  equations 


F  -  c  P2/(P1  cos  G)  (6a) 

G  a  c  P3/(PX  cos  F)  (6b) 


where  c  is  a  constant  and  F and  G  are  the  complements  of  the 
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directional  angles.  Accordingly, 

w  =  V  sin  F  (7a) 

v  =  V  sin  G  (7b) 

u  »  V(cos  F) (cos  G)  (7c) 

The  components  of  the  wind  vector  are  described  more 
closely  by  equations  (3) ,  (6) ,  and  (7)  than  by  equations 
(3) ,  (4)  ,  and  (5) .  When  using  (6) ,  an  iterative  procedure 
is  used  to  solve  for  F  and  G  which  are  then  used  in  equa¬ 
tions  (7)  . 

b.  Pressure  transducers 

Capacitive  pressure  transducers  manufactured  by 
2/ 

Datametncs,  Inc.-'  were  chosen  for  the  pressure  measure¬ 
ment.  The  gains  of  the  signal  conditioners  can  be  selected 
to  provide  full  scale  outputs  (+  5.0V)  for  differential 

pressures  of  10,  20,  30,  GO,  100,  200,  300,  600,  1000, 

_2 

2000,  3000,  6000,  10,000  dynes  cm  .  The  transducer  has 

a  maximum  nonlinearity  of  about  +  0.1%,  zero  drift  of 

10~5  of  maximum  range  per  degree  Celsius  and  sensitivity 

-2 

change  of  2  x  10  %/C. 

c.  Frequency  response 

The  frequency  response  and  phase  shi::t  of  a 
pressure  transducer  connected  by  tubing  to  a  fluctuating 
pressure  has  been  described  by  Iberall  (1950) ,  whose 


U  Datametrics  Incorporated,  87  Beaver  Street,  Waltham, 
Mass.  (Model  511-8  Barocel) . 
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analysis  was  basic  to  our  system  design.  The  response 

of  the  transducer  is  controlled  by  the  size  and  length  of 

the  tubing  and  the  effective  internal  volume  of  the 

transducers.  The  transducers  used  in  our  system  were  a 

special  design  which  used  a  stif fer- than -normal  diaphragm 

3 

and  a  reduced  internal  volume  of  1.6  cm  to  improve  the 
frequency  response  of  the  system.  The  spherical  probe 
was  connected  by  approximately  43  cm  of  approximately 
1 . 5 -mm  I.D.  tubing  to  the  pressure  transducer;  tests 
showed  this  tubulation  optimized  the  system  performance. 

The  frequency  response  and  phase  shift  of  the  system 
were  checked  by  producing  known  sinusoidal  pressure  dif¬ 
ferences  at  various  frequences  between  appropriate  ports 
on  the  surface  of  the  pressure  spheres  and  monitoring 
the  amplitude  and  phase  of  the  transducer  output.  In¬ 
stead  of  attempting  to  produce  the  pressure  differences 
between  ports  on  a  single  sphere,  two  identical  spheres 
were  placed  in  separate  pressure  chambers  with  tubing 
connecting  appropriate  ports  to  the  pressure  transducers. 
Equal  pressure  fluctuations,  180  degrees  out  of  phase 
with  each  other,  were  produced  in  the  two  chambers  by 
pistons  which  were  closely  coupled  to  the  chambers. 

The  pistons  were  dri'  en  by  a  variable  speed  motor  and 
the  phase  of  the  p:  assure  fluctuation  was  determined  by 
optically  sensing  the  position  of  the  Scotch  yoke  piston 
drive.  Typical  amplitude  and  phase  shift  characteristics 
of  the  3-cm  anemoclinometer  and  pressure  transducer  are 
shown  in  Table  1.  The  response  was  limited  by  the  tubing 
used  to  construct  the  anemoclinometers  and  could  be 
improved  by  redesigning  the  pressure  sphere,  tubulation, 
and  transducer  system  for  optimum  performance. 
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Table  1.  Anemometer  system  frequency  response  and  phase 
shift. 


frequency 

Hz 

relative 

amplitude 

phase  shift 
(degrees) 

1 

1.00 

0 

5 

1.02 

18 

10 

1.05 

48 

15 

1.00 

76 

20 

0.84 

100 

25 

0.63 

126 

30 

0.47 

145 

d.  Field  installation 

For  field  measurements,  the  pres sure- sphere  anemometer 
is  mounted  on  a  2.5-cm  diameter  mast  (Fig.  3)  at  the  desired 
height.  The  anemometer  is  oriented  with  the  shaft  axis 
parallel  to  the  anticipated  direction  of  mean  flow.  The 
three  pressure  transducers  are  housed  in  a  temperature- 
controlled  (_+  0.2C)  box  which  is  an  integral  part  of  the 
mounting  assembly  located  at  the  opposite  side  of  mast  to 
the  pressure  sphere.  The  temperature  control  provides 
the  required  zero  stability.  The  whole  assembly  can  be 
moved  co  different  levels  on  the  mast  or  completely 
removed  as  one  unit  without  disconnecting  the  pressure 
transducers  from  the  pressure  sphere.  The  pressure 
transducers  are  connected  by  150  m  of  cable  to  their 
power  supply  and  signal  conditioners  which  are  housed  in 
a  2.5  x  6  m  air-conditioned  instrument  trailer. 

The  masts  are  on  pivot  points  and  supported  higher 
up  by  guy  wires  attached  to  bearings.  This  arrangement 


Fig.  3.  Anemometer  assembly  on  its  mast. 
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allows  the  mast  to  be  rotated  so  that  the  probe  can 
be  oriented  easily  into  the  mean  wind.  When  the  data 
discussed  below  were  obtained,  the  anemometer  was 
rotated  in  azimuth  manually  into  the  mean  wind;  at  the 
beginning  of  each  half-hour  run  the  orientation  was 
adjusted  to  the  position  of  the  mean  wind  for  the  previous 
half  hour.  Since  then,  a  motor  assembly,  controlled  by 
the  =  bpuv  output  of  the  wind  probe,  has  been  used  to 
continuously  but  slowly  adjust  the  position  of  the  probe 
into  the  wind.  The  orientation  of  the  mast  is  monitored 
through  the  output  of  a  potentiometer  attached  to  the 
base  of  the  mast  and  is  included  in  the  calculation  of 
the  components  of  the  wind  vector. 

3.  Data  handling 

Since  the  sensors  respond  to  frequencies  as  high  as 
30  Hz,  a  large  amount  of  data  must  be  analyzed  if  the 
system  is  operated  over  extended  periods  of  time.  Storing 
large  quantities  of  data  under  field  conditions  is  costly 
and  often  results  in  a  serious  reduction  in  data  quality. 

In  addition  it  is  highly  desirable  that  some  data  analysis 
capability  be  available  at  the  experimental  «ite  so  that 
the  experiment  can  run  efficiently  and  urumentation 
faults  detected  as  scon  as  possible.  After  a  careful  study 
of  the  available  alternatives  we  elected  to  drastically 
reduce,  by  digital  on-line  computation,  the  quantity  of 
data  to  be  stored  to  the  point  where  it  could  be  typed 
out  in  table  form  by  a  typewriter  or  stored  on  paper  tape. 
In  1967,  this  amounted  to  a  data  reduction  of  approxi¬ 
mately  18000:1. 

In  1957  the  data  analysis  was  performed  on  an  EMR  6020 
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computer-  and  later  on  the  faster  and  smaller  EMR  6130, 

The  1967  system  included  a  Raytheon  A-D  converter,  6020 
computer  and  model  33  teletype  with  paper  tape  reader  with 
punch.  The  6130  system  includes  an  EMR  2701  converter, 
and  a  higher  speed  paper  tape  reader  and  punch  in  aduition 
to  the  model  33  teletype. 

Five  channels  of  analogue  data  were  obtained  at  each 
of  three  sites  to  give  a  total  of  15  channels.  At  each 
site  three  channels  represented  the  three  pressure  differ¬ 
ences  P^,  P2#  and  the  other  two  channels  represented 
a  fast  response  resistance  thermometer  (Wesely,  et  al., 
1969)  and  a  fast  response  barium  fluoride  relative  humid¬ 
ity  element  (Jones,  1967) .  The  velocity  components  of 
the  wind  vector  were  calculated  using  equations  (6)  and 
(7)  and  means,  squares,  and  crossproducts  of  the  five 
parameters  (u,  v,  w,  T,  e)  were  calculated,  where  T,  and 
e  were  the  temperature  and  vapor  pressure  respectively. 

The  complete  operation  (i.e.  15  channels  of  analogue 
to  digital  conversion  and  the  data  analysis)  was  repeated 
40  times  per  second.  At  the  end  of  each  half-hour  sampling 
period  the  necessary  scaling  operations  were  performed  and 
the  outputs  were  fceletyped.  Approximately  2.5  minutes  of 
each  half-hour  period  were  required  for  output  and  no  data 
were  collected  during  this  time. 

This  data  system  has  proven  to  be  a  very  efficient 
and  powerful  research  tool  and  it  is  felt  that  the  success 
achieved  with  the  anemometer  system  would  not  have  been 
possible  if,  alternatively,  data  storage  equipment  had 
been  selected. 


3/ 

—  Electro-mechanical  Research,  Inc.,  0001  Bloomington 
Freeway,  Minneapolis,  Minn. 
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4 .  Tests  of  anemometer  system 

A  complete  description  of  the  experimental  area  may 
be  found  in  Brooks  and  Pruitt  (1966) .  A  plan  of  the 
field  site  is  given  in  Fig.  4,  showing  the  heights  and 
spatial  arrangement  of  our  three  anemometers  with  respect 
to  the  6-meter ,  Davis  shear  stress  lysimeter  and  the 
triangular  array  of  masts  installed  by  Dr.  C.  R.  Stearns, 
University  of  Wisconsin  Department  of  Meteorology.  These 
masts  carried  cup  anemometer  and  aspirated  dry-  and  wet- 
bulb  thermometers.  The  surface  was  uniform  aita  fescue, 

5  to  10  cm  high,  which  was  periodically  mown. 

Wind  velocity  measurements  with  our  pressure  probe 
are  compared  with  cup  anemometer  data  and  our  eddy- 
correlation,  shear  stress  measurements  are  compared  with 
both  shear  stress  lysimeter  data  and  shear  stresses 
obtained  by  Dr.  Stearns*  preliminary  analysis  of  his 
vertical  profiles  of  wind  velocity  and  of  temperature 
(KEYPS-type,  diabatic  profile  analysis) .  In  addition, 
a  graphical  description  of  the  vertical  fluctuations 
of  wind  velocity  as  a  function  or  the  stability  param¬ 
eter  z/L  is  presented. 

a .  Comparison  of  wind  velocity  measurements 

The  on-line  computer  program  which  was  used  to 
analyze  our  anemometer  data  included  the  calculation  of 
the  horizontal  wind. 

=  (u2  +  v2)*5 

where  u  and  v  are  the  instantaneous  values  cf  the  hori¬ 
zontal  components  of  the  wind  vector.  The  value  of  vT 
is  primarily  dependent  upon  P^,  as  given  in  (2) ,  and 
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since  u  generally  is  much  larger  than  either  v  or  w, 
errors  associated  with  '.he  measurement  of  wind  angles 
calculated  from  (6)  do  not  seriously  degrade  the 
estimate  of  V,,.  The  good  agreement  between  V~  and  cup 
anemometer  measurements  presented  in  Fig.  5,  demonstrate 
the  accuracy  of  pressure-sphere  measurements  of  u. 

b .  Comparison  of  shear  stress  measurements 

Shear  stress  measurements  obtained  with  che  pressure- 
sphere  anemometer  are  compared  with  data  from  the  shear- 
stress  lysimeter  and  from  analysis  of  the  wind  profiles. 

The  data  obtained  on  May  2,  2,  4,  and  5  are  presented  in 
Figs.  6  and  7.  The  pressure-sphere  anemometer  data 
represent  the  average  of  measurements  available  at,  the 
three  sites.  The  shear  stress  data  from  the  three  wind 
profile  sites  also  were  averaged.  The  Davis  shear  stress 
lysimeter  independently  measures  the  north-south,  and 
east-west  components  of  the  surface  shear  stress  and  the 
data  used  were  computed  by  W.  0,  Pruitt  as  the  vector  sum 
of  the  half-hour  means  of  these  components. 

Agreement  among  the  three  methods  is  satisfactory  even 
though  the  aerodynamic  analysis  generally  provides  some¬ 
what  larger  values  than  the  other  two  methods.  This  dis¬ 
crepancy  appears  unduly  large  on  May  4  and  5.  The  average 
zq  value  computed  from  the  wind  profiles  is  0.95  cm  for 
May  2  and  3,  and  about  1.4  cm  for  May  4  and  5.  For  the 
latter  two  days  new  estimates  of  the  shear  stress  were 
calculated  via  a  drag  coefficient  using  =0.95  cm, 
a  KEYPS  diabatic  correction  and  the  cup  anemometer  wind 
velocity  at  80  cm.  The  results  of  this  calculation  are 
more  consistent  with  the  comparisons  on  May  2  and  3. 
Calculations  indicate  that  best  agreement  between  drag 
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Fig.  6.  Comparison  of  shear  stress  measurements  on 
May  2  and  3,  1967. 
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Fig.  7.  Comparison  of  shear  stress  measurements  on 
May  4  and  5 t  1967. 
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coefficient  and  eddy  correlation  determinations  would 
have  been  obtained  by  using  zq  **  0.7  cm. 

c .  Standard  deviation  of  the  vertical  wind 

Ratios  of  the  standard  deviation  of  the  vertical 
component  of  wind  velocity  to  the  friction  velocity,  u*, 
as  measured  at  the  two  1-m  sites  and  one  4-m  site,  are 
plotted  in  Fig.  8.  The  comparison  of  horizontal  wind 
measured  with  our  anemometer  system  and  with  cup  anemom¬ 
eters  indicates  that  our  pressure-sphere  anemometer 
measures  the  u-component  of  wind  velocity  accurately. 
Accordingly  the  ratio  r?w/u*  would  vary  as  the  square  root 
of  a  constant  percentage  error  .in  the  measurement  of  w. 
This  is  not  a  very  sensitive  test  of  the  measurement 
of  the  vertical  component  of  velocity  since  the  error 
in  w  would  be  twice  that  in  crw/u*,  but  our  value  of 
1.25  for  *?w/utt  under  neutral  conditions  is  the  same  as 
that  derived  by  Panofsky,  et  al.,  (1967)  and  close  to 
the  value  of  1.33  predicted  by  Lettau  (1968) .  Over  one 
hundred  additional  data  points  were  collected  over  snap 
bear  s  (zq  -  4  cm)  in  1968  at  Hancock,  Wisconsin.  The 
data,  averaged  over  stability  ranges  of  z/L  =  -  0.^5  to 
0.25,  are  also  presented  in  Fig.  8  and  are  very  similar 
to  the  Davi3  data.  For  ready  comparison  with  the 
panofsky  (1967)  and  panofsky  et  al.  (1967)  the  curve 
[1  -  ( z/L) /s] ^  is  plotted. 
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ABSTRACT 

A  three-dimensional  pressure- sphere  anemometer  and 
fast  thermometer  system  (P.S.A.T.)  was  used  to  measure 
vertical  heat  flux  density  in  the  atmospheric  surface 
layer  at  one  to  four  meters  above  alta  fescue  and  snap 
beans.  Good  agreement  with  independent  measurements 
was  obtained,  which  shows  that  the  P.S.A.T.  is  suf¬ 
ficiently  small  and  has  adequately  high  frequency 
response  and  accuracy  for  eddy-correlation  measurements 
within  one  meter  of  the  surface.  Also  obtained  with 
the  P.S.A.T.  were  (u't' ) /{w't' ) ,  ru>T,  rw>T»  and  oT/T* 
and  their  dependence  upon  stability.  When  the  atmos¬ 
phere  was  thermally  stable,  slow  wave  motions  frequently 
increased  <?T  even  though  turbulent  mixing  was  lacking. 
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1 .  Introduction 

The  turbulent  vertical  heat  flux,  H,  in  the  atmos¬ 
pheric  surface  layer  ever  a  horizontally  uniform  surface 
can  be  determined  from 


j  H  =  p  w't'  (1) 

i 

where  o  is  the  air  density,  c^  is  the  specific  heiat  of 
air,  w  is  the  vertical  wind  velocity,  and  T  is  the  air 
temperature.  The  bar  denotes  a  time  average  and  the  prime 
denotes  an  instantaneous  deviation  from  the  time-averaged 
quantity.  The  major  difficulty  with  making  eddy  correlation 
measurements  of  turbulent  heat  transport  is  in  measuring 
the  vertical  wind.  This  requires  an  accurate  and  stable 
anemometer  that  measures  the  wind  components  with  a  suf¬ 
ficiently  high  frequency  response  for  use  close  to  the 
surface  where  the  fetch  requirements  are  minimum.  At 
present,  the  most  promising  anemometers  are  sonic  anemom¬ 
eters,  either  pulsed  wave  (Mitsuta,  1966)  or  continuous 
wave  (Kaimal,  et  al.,  1968),  and  the  pressure-sphere 
anemometer  (Thurtell,  et  al.,  1969).  The  pressure-sphere 
anemometer  is  smaller  than  sonic  anemometers  and  thus  can 
|  be  used  closer  to  the  surface  Where  eddies  are  smaller. 

This  paper  describes  the  measurement  of  turbulent 
heat  transport  with  the  pressure-sphere  anemometer  and  a 
small,  fast- response,  resistance  thermometer.  Measure¬ 
ments  of  heat  flux  above  alta  fescue  are  compared  with 
independent  measurements  made  by  others  at  the  University 
of  California  at  Davis  as  part  of  the  1967  Cooperative 
!  Field  Experiment  sponsored  by  the  Atmospheric  Sciences 

j  Laboratory,  U*  S.  Army  Electronics  Command,  Ft.  Huachucs, 

i  Arizona.  Also  presented  arc  measurements  of  heat  flux 

I 
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made  in  1968  over  snap  beans  at  the  University  of  Wisconsin 
Hancock  Experiment  Farm.  A  summary  of  the  standard  devi¬ 
ation  of  temperature,  c/^,  divided  by  the  dimensionless 
temperature  scale,  T*,  and  of  correlation  coefficients 
for  wind  and  temperature  also  is  given  for  measurements 
over  the  above  two  surfaces. 

2.  Equipment,  sites,  and  comparison  measurements. 

a  •  Pressure  sphere  anemometer-thermometer  assembly 

A  fine-wire-resistance  thermometer  (Wesely,  et  al.# 
1969)  was  mounted  parallel  to  the  horizontal  ports  of  the 
pressure  sphere  (Thurtell.  et  ai.#  1969),  as  shown  for 
a  3-cm  diameter  sphere  in  Fig.  1.  The  closest  edge  of 
the  thermometer  was  about  1.25  cm  from  the  3-cm  anemo- 
clinometer  and  about  2.5  cm  from  the  8-cm  anemoc] inometer . 
The  thermometer  was  placed  so  that  the  sensitive  element 
was  slightly  upwind  of  the  leading  edge  of  the  pressure 
sphere;  tests  showed  that  this  forward  placement  was 
necessary  to  prevent  thermal  modification  of  the  air  that 
flowed  to  the  thermometer  past  the  large  thermal  mass  of 
the  sphere.  The  thermometer  was  outside  the  angle  of 
acceptance  of  the  anemoclinometer  and  tests  showed  that 
the  flow  patterns  around  the  sphere  were  not  significantly 
affected . 

When  the  1967  data  were  obtained,  the  anemometer  was 
rotated  in  azimuth  manually  in.o  the  moan  wind  at  the 
beginning  of  each  half-hour  run,  During  J.968,  a  motor 
assembly  controlled  by  the  anemometer,  rotated  the  mast 
to  point  tne  anemometer  into  t .ho  wind.  The  azimuth 
rotation  of  the  mast  v:as  monitored  with  a  potentiometer 
attached  to  the  base  of  the  mast  and  was  included  in  the 


Fig.  1.  Arrangement  of  thermometer  with 
the  anemocl.inometer . 


calculation  of  the  components  of  the  wind  vector. 

b .  Data  handling 

The  current  through  the  thermometer  was  kept  nearly 
constant  at  0.3  ma  by  its  bridge,  which  was  located  about 
5  m  from  the  thermometer.  The  bridge  output  was  fed 
directly  into  a  floating  differential  amplifier  with  a 
1000  gain,  to  provide  a  signal  with  a  temperature  sensi¬ 
tivity  of  about  0.6C  v”1.  The  amplifier  output  was 
transmitted  through  150  m  of  cable  to  the  instrument 
trailer.  The  thermometer  and  anemometer  signals  were 
fed  to  a  scanner-converter  and  an  EMR  computer  as  de¬ 
scribed  by  Thurtel.l,  et  al.f  (1969).  The  sampling  rate 
was  40  sec  1  in  1967  and  150  sec”1  in  1966. 

The  outputs  of  the  thermometer  bridges  were  filtered 
in  the  amplifiers  to  match  the  phase  shifts  and  response 
of  the  pressure- sphere  anemometer  and  a:  so  to  avoid  high 
frequency  noise.  The  response  of  the  two  systems  are 
shown  in  Fig.  2.  The  curves  for  the  anernoclinometers 
are  roughly  representative  of  the  vertical  wind  component. 

c .  Site  description 

A  description  of  the  site  of  the  1967  Cooperative 
Field  Experiment,  and  our  instrument  locations  as  well  as 
the  locations  of  other  relevant  instruments  may  be  found 
in  Thurtell,  et  al.,  (1969). 

The  1968  measurements  at  Ilancoc!;,  Wisconsin  were  on 
a  100  x  160  m  field  of  snap  beans  planted  in  rows  spaced 
at  90  cm.  The  snap  beans  were  about  30  cm  high  and 
provided  about  50%  cover  over  Plainfield  sand.  The 
fetch  was  60  m  to  the  north,  50  m  to  the  east  and  west 
and  100  m  to  the  south.  Beyond  these  boundaries  to  the 
south  was  alfalfa  extending  for  150  m  to  a  15-m  high 
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Fig.  2.  Frequency  response  and  phase  shift  of  the 
thermometer  system  and  of  the  anemometer  system 
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woods  and  to  the  west  was  an  alfalfa  field  extending 
100  m  to  a  10-m  high  shelter  belt,  to  the  northwest  were 
low  crops  extending  200  m  to  a  shelter  belt,  and  to  the 
east  was  alfalfa  extending  300  m  to  a  woods.  The  wind 
was  predominate.! y  from  the  south  and  west  during  the 
tests. 

d.  Independent  measurements  of  sensible  heat  flux 

During  the  1967  Cooperative  experiment,  Dr„  C,  R. 
Stearns,  University  of  Wisconsin,  measured  wind,  tempera¬ 
ture,  and  vapor  pressure  profiles  at  three  locations  in 
a  triangular  array.  At  the  same  locations,  he  measured 
net  radiation  and  soil  heat  flux  density  for  energy 
balance  calculations.  The  sensible  heat  flux  was  calcu¬ 
lated  from  the  energ  balance  using  Bowen1 s  ratio, 

Q  =  y&T /AeP  determine  1  from  vertical  temperature  and 
vapor  pressure  differences  measured  over  the  same  height 
intervals  within  120  cm  of  the  surface.  An  aerodynamic 
calculation  of  the  sensible  heat  flux  also  was  made  using 
the  wind  and  temperature  profiles  to  find  the  shear 
stress  with  a  KEYPS-type  analysis  and  then  using  simi¬ 
larity  (Ku  =  K.,)  and  the  profiles  to  find  the  heat  flux. 
Dr.  Stearns  supplied  us  data  from  both  analyses. 

The  University  of  Calif ornia-Davis  group  measured 
the  evaporation  with  a  o  -m  diameter  weighing  lysimetcr 
(Pruitt  and  Angus,  1960) .  In  addition,  they  measured 
net  radiation  and  soil  heat  flux  near  the  lysimeter. 

The  sensible  heat  flux  density  was  calculated  by  dif¬ 
ferencing  the  energy  balance  terms  as  II  =  K^-G-E,  The 
Davis  group  also  measured  the  sensible  heat  flux  directly 
with  an  Evapotron  (Dyer  and  Mahei ,  1965).  Both  measure¬ 
ments  were  supplied  to  us  by  Dr.  W.  0.  Pruitt,  University 
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of  California-Davis . 

The  University  of  Washington  group  measured  the 
sensible  heat  flux  both  with  a  one-dimensional  sonic 
anemometer  thermometer  (Kaimal  and  Businger,  1963)  and 
with  a  three-dimensional  unit  (Mitsuta,  1966) .  These 
data  were  supplied  by  Dr.  J.  A.  Businger. 

The  comparison  data  at  Hancock,  Wisconsin  were 
obtained  by  differencing  measurements  of  net  radiation, 
soil  heat  flux  density  and  evaporation.  The  evaporation 
was  measured  with  a  2.1  x  5.5  m  weighing  lysimeter 
(Black  et  al.,  1968),  the  net  radiation  was  measured 
with  a  large  Funk  radiometer,  and  the  soil  heat  flux 
was  measured  with  soil  heat  flux  plates  (Fuchs  and  Tanner, 
1968)  and  integrating  soil  thermometers  (Tanner,  1958) . 

3.  Heat  flux  density  comparisons 

During  the  1967  Cooperative  Field  Experiment,  fetches 
were  easily  in  excess  of  100  m,  except  for  small  changes 
in  elevation,  since  the  wind  was  predominately  from  the 
south  and  southwest  where  fields  had  similar  vegetation 
and  roughness. 

Heat  flux  estimates  by  the  pressure-sphere  anemometer 
and  thermometer  system  (called  the  P.S.A.T.  hereafter)  were 
averaged  from  two  1-m  high  units  and  one  4-m  high  unit  to 
cive  the  results  shown  in  Fig.  3.  There  was  no  systematic 
difference  of  heat  flux  measured  at  the  two  heights  except 
from  1615  to  2015  on  April  27,  when  the  data  from  the 
higher  nast  were  discarded.  Heat  flux  data  from  a  three- 
dimensional,  sonic  anemometer-thermometer  at  four  meters 
above  the  surface  and  one-dimensional,  sonic  anemometer- 
thermometer  2.2  m  high  agreed  well  with  the  P.S.A.T.; 
the  scatter  of  estimates  at  our  three  different  P.S.A.T. 
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Comparisons  of  heat  flux  density  estimates  at 
Davis,  California  for  April  27,  May  2,  May  3, 
and  May  4  and  5  of  1967. 
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sites  frequently  is  of  the  same  order  as  the  difference 
between  our  data  and  that  of  the  sonic  anemometer- 
thermometer  . 

The  eddy-flux  from  the  Evapotron  is  shown  on  Pig,  3 
for  May  4  and  5.  The  wide  fluctuations  may  have  been  due 
to  the  averaging  process  to  remove  the  me^n  wind  and 
temperature  terms  since  a  time  constant  of  only  one 
minute  is  used  in  this  system. 

Several  indirect  estimates  of  sensible  heat  flux  are 
also  shown  in  Fig.  3.  The  energy  balance  estimates 
obtained  from  differencing  the  energy  balance  H  ■  Rn~G-E 
appear  high  during  the  day  and  low  at  night.  Since  |  Rn( 
and  | Ej  are  much  larger  than  |h|c  a  small  relative  error 
in  these  terms  could  produce  a  large  relative  error  in 
!  H|  . 

The  results  from  the  Bowen  ratio  energy  balance  and 
those  from  the  aerodynamic  method  are  the  averages  of  heat 
flux  data  from  three  sites.  These  methods  are  nearly  inde¬ 
pendent,  but  not  completely  so,  because  they  use  the 
same  temperature  profiles.  Both  methods  show  remarkably 
good  agreement  with  the  P.S.A.T. 

In  Fig.  4  is  shown  the  average  of  heat  flux  estimates 
at  two  P.S.A.T.  sites.  Both  sites  were  117  cm  above  the 
soil  surface  until  1030  when  one  site  was  moved  to  210  cm 
above  the  surface.  Since  estimates  of  heat  flux  by  the 
P.S.A.T.  at  210  cm  from  the  soil  surface  were  not  system¬ 
atically  different  from  the  117-cm  high  site,  fetch  was 
considered  adequate.  On  another  day,  we  compared  measure¬ 
ments  with  one  P.S.A.T.  at  75  cm  and  the  other  at  117  cm 
and  found  no  systematic  differences. 

The  energy  balance  estimate  of  heat  flux  leads  the 
P.S.A.T.  estimate  in  the  morning.  This  was  probably  due 
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to  a  time  lag  in  the  evaporative  flux  caused  by  unrepre¬ 
sentative  heat  storage  in  the  lysimeter  (Black,  et  al., 

1968) .  This  also  could  have  caused  an  overestimate  of 
the  magnitude  of  the  heat  flux  after  sunset.  The  low 
value  ^  i.  1015  was  caused  by  an  unexplainably  large 
estimate  of  evaporative  flux. 

4.  Temperature  structure 

When  the  data  used  in  this  section  were  collected, 

V.  and  t  were  constant  with  height,  within  the  accuracy 

of  our  measurements;  thus  we  c  m  use  H  and  t  as  scaling 

factors  a. 3  described  by  Monin  and  Obukhov  (1954) .  They 

define  a  dimensionless  height  ratio  z/L  where  z  is  the 

3 

height  from  the  surface  and  L  =  -u*  0c^T/(kgH)  [u*  =  (t/q)  • 

is  the  friction  velocity;  k  =:  0.428  is  Kerman' s  constant, 

-2 

and  g  =  980  cm  sec  ] .  The  relationships  obtained  between 
our  measurements  cf  z/L  and  our  measurements  of  the  corre¬ 
lation  coefficients  T  and  ^  ,  and  of  the  ratio 
(u't' ) /(w't' )  are  given  in  Fig.  5.  Fig.  7  shows  the 
relation  of  z/L  to  a  dimensionless  standard  deviation  of 
temperature,  a^T*  [aT  is  the  standard  deviation  of  air 
temperature  and  T*  =  -K/ (ko c^u^) ] . 

It  appears  that  (u't* ) / (w't' )  «  4  for  z/L  =  0.1  and 
«  2.5  for  z/L  =  -0.05.  The  large  scatter  indicates  that 
more  meaningful  results  might  have  been  obtained  from 
sampling  periods  shorter  than  the  30  min  used.  For 
instance,  Zubkov skii  and  Tsvang  (1966)  obtain  less 
scatter  by  using  running  means  of  the  winds  and  tempera¬ 
tures  from  electrica]  filters  with  time  constants  of 
100  sec  and  80  sec,  respectively.  Fig.  5  shows  that  air 
temperatures  are  more  closely  coupled  with  horizontal 

winds  than  with  vertical  w inds  since  |r  ,J>|r  . 
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Fig.  5.  Correlation  coefficients  of  wind  and  tempera¬ 
ture  and  (u'T* ) /(wT7 )  as  a  function  of  stabil 
ity. 
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This  is  especially  true  for  stable  conditions.  As  shown 
in  Fig.  6,  the  fluctuations  of  air  temperature  and 
vertical  wind  during  unstable  conditions  are  much  larger 
and  faster  than  during  stable  conditions;  however,  it 
has  been  observed  that  slow  wave  motion  (not  evident  in 
Fig.  6)  frequently  occurs  at  night  when  the  wind  speed 
is  low.  Then  the  temperature  at  a  stationary  height 
in  a  highly  stratified  atmosphere  fluctuates  as  much  as 
3C  every  50  to  200  sec,  the  period  of  the  slow  waves. 

This  oscillation  substantially  increases  aT  ard  probably 
accounts  for  some  small  values  of  r  _  and  |r  |  for  large 
positive  values  of  z/L. 

In  Fig.  7,  a,p/T*  is  plotted  against  z/L  and  aprears 
to  scale  well  for  unstable  conditions,  except  near  z/L  =  0, 
v/here  T*  =  0.  A  function  suggested  by  Dyer  (1965)  is 
drawn  following  Panofsky,  et  al.,  (1967).  Data  from 
Russian  sonic  anemometers  and  resistance  thermometers 
(Mordukhovich  and  Tsvang,  1966)  and  data  from  a  one¬ 
dimensional  sonic  a nemometer- thermometer  'Businger  et  al., 
1967)  are  also  included  in  Fig.  7.  The  P.3.A.T.  data 
agrees  well  with  the  Russian  data,  but  appears  lower 
than  the  data  summarized  by  Panofsky  et  al.,  (1967). 

The  large  scatter  for  stable  conditions  may  be 
caused  in  part  by  small  absolute  errors  in  t  and  H,  since 
both  t  and  H  are  about  ten  times  smaller  at  night  than 
during  the  day;  however,  slow  wave  motion  may  increase 
t  without  increasing  the  heat  flux  enough  to  keep 
rrf|1/T*  from  increasing  whenever  these  large-scale  dis¬ 
turbances  occur.  Since  Mordukhovich  and  Tsvang  (1966) 
use  a  running  mean  of  temperature  with  a  time  constant 
of  £0  sec,  temperature  oscillations  with  peri»xis  longer 
than  20  sec  are  substantially  attenuated,  causing  their 
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Fig.  6.  Fluctuations  of  T  and  nuw  for  stable  condition 
at  2.0  m  (bottom)  above  bare  Plainfield  sand 
at  Ha ncock ,  V? i 3 . 
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Fig.  7.  Standard  deviation  of  dimensionless  tempera 
ture  as  a  function  of  stability. 


42 


estimates  of  r^/T*  to  have  less  scatter  and  be  lower 
than  our  estimates.  Measurements  during  stable  con¬ 
ditions,  when  the  wind  speeds  were  at  least  2.5  m  sec”'1' 
at  1  m,  have  less  scatter:  mixing  is  probably  adequate 
then  to  prevent  domination  by  large-scale  disturbances. 
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ABSTRACT 

Eddy  correlation  measurements  of  water  vapor  flux 
density  have  been  made  using  a  barium  fluoride  film 
humidity  sensor.  During  morning  and  evening  periods, 
good  agreement  was  obtained  between  eddy  correlation 
data  and  two  independent  methods.  Serious  disagreement 
between  measurements  occurred  only  when  the  humidity 
sensor  was  operating  within  a  poorly  defined  portion  of 
the  calibration  curve  which  was  not  suited  to  on-line 
calculations.  The  results  indicate  that  the  humidity 
sensor  could  be  modified  to  allow  operation  at  all 
times  within  well  defined  segments  of  the  calibration 
curve  and  permit  successful  eddy  correlation  vapor  flux 
measurements  within  one  meter  of  the  surface.  (Key  Vlords 
Humidity  sensor;  eddy  correlation;  vapor  flux) 
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Introduction 

Of  the  micrometeor ological  methods  currently  avail¬ 
able  for  determining  evaporation,  the  eddy  correlation 
approach  is  most  satisfying  since  it  requires  the  least 
number  of  basic  assumptions.  The  equation  which  describes 
the  evaporation  as  latent  heat  flux  density,  may  be 
written  as 


Fv  =  X  [q  w  +  q'w']  (1) 

where  >.  is  the  .latent  heat  of  vaporization,  q  is  vapor 
concentration  (absolute  humidity)  and  w  is  the  vertical 
wind  velocity.  The  overbars  indicate  time  averages  and 
the  primes  indicate  fluctuations  about  the  mean.  The 
surface  evaporation,  E,  will  be  equal  to 

E  =  >.  [q7^]  (2) 

when  w  is  equal  to  zero. 

Although  eddy  correlation  measurements  of  sensible 
heat  flux  have  beer,  made  [Kaimal  and  Businger,  1963? 
Businger,  et  ai.,  1967?  Wesely,  Thurtell,  and  Tanner, 
1969]  evaporation  measurements  have  been  limited  by 
slow  humidity  sensors.  Dyer  and  hicks  [1967]  and 
Goddard  and  Pruitt  [1966],  using  fine-wire  psychrometers, 
found  that  measurement  at  four  meters  was  necessary 
where  larger  and  slower  eddies  could  be  recorded  by 
these  relatively  slow  elements.  At  these  elevations, 
however,  storage  and  advection  errors  occurred 
unless  there  was  a  very  long  fetch.  In  order  tc  work 
closer  to  the  ground,  we  have  investigated  the  possible 
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use  of  a  rapid-response  barium  fluoride  film  humidity 
sensor  [Jones,  1967],  Bean  and  Florey  [1968]  report  on 
the  use  of  this  sensor  for  measuring  evaporation  at  two 
meters  above  Lake  Hefner.  However,  their  system  was 
limited  by  the  relatively  slow  response  of  an  anemomete.r- 
bivane  and  not  by  the  humidity  element's  response.  We 
believe  that  in  association  with  a  fast  response  wind 
vector  sensor  the  barium  fluoride  film  humidity  sensor 
can  allow  measurement  of  evaporation  considerably  closer 
to  the  surface. 


Instrumentation  and  Methods 

Barium  fluoride  film  humidity  sensor .  The  barium 
fluoride  humidity  sensor  consists  of  a  glass  plate  of 
approximately  10  x  2  x  0.16  cm  on  which  a  0.3vi -thick 
film  of  barium  fluoride  has  been  evaporated  over  closely- 
spaced,  evaporated  chromium  electrodes.  The  electrical 
resistance  is  measured  between  the  electrodes.  Jones 
[1967]  reports  in  detail  on  these  sensors  and  their 
properties. 

The  particular  elements  used  in  the  present  work  are 
calibrated  by  determining  their  resistances  over  a  series 
of  known  relative  humidity  solutions  from  12  to  97%. 

Plots  of  the  logarithm  of  sensor  resistance  against  rel¬ 
ative  humidity  consist  of  three  linear  segments.  Un¬ 
fortunately,  sensor  calibration  is  not  stable  for  an 
unlimited  time  and  degrades  substantially  over  a  period 
of  several  months.  Calibration  curves  indicating  changes 
over  time  are  shown  in  Figure  1. 

In  order  to  use  relative  humidity  measurements  in 
the  eddy  correlation  method,  a  reference  temperature  must 
be  measured.  During  preliminary  tests  at  Davis,  California 
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Calibration  curves  for  a  barium  fluoride 
film  humidity  sensor  showing  changes  with 
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in  1967,  we  used  the  air  temperature;  the  results  showed 
that  the  sensor  film  temperature  should  be  used  as  the 
reference  for  converting  relative  humidity  to  absolute 
humidity  and  must  be  monitored.  Accordingly,  a  127y. 
micro-bead  thermistor  is  cemented  with  a  very  small 
amount  of  clear  epoxy  to  the  center  of  the  sensitized 
surface.  A  linearized  bridge  is  used  with  the  thermistor 
for  the  surface  temperature  measurement. 

The  sensors  also  have  electrical  leads  cemented  to 
them.  Vapors  from  the  epoxy  used  to  attach  the  leads  and 
thermistors  caused  an  immediate  calibration  shift,  and 
it  is  possible  that  the  drifts  shown  in  Figure  1  were 
accelerated  by  the  early  exposure  to  organic  vapors. 

Electronic  circuitry.  A  block  diagram  of  the  cir¬ 
cuit  used  with  the  barium  fluoride  humidity  sensor  and  its 
associated  thermistor  is  shown  in  Figure  2.  The  two  most 
important  features  of  the  system  are  the  logarithmic  ampli¬ 
fier  and  the  phase  adjustment.  The  logarithmic  amplifier 
provides  an  output  voltage  that  is  linear  with  relative 
humidity  as  shown  in  Figure  1.  The  phase  adjustment  is 
necessary  because  at  low  humidities,  and  with  very  high 
sensor  resistances,  there  is  significant  capacitive  react¬ 
ance.  The  phase  is  adjusted  to  null  the  capacitive  react¬ 
ance  while  the  sensor  is  at  very  low  humidity  over  a 
desiccant  ■  no  further  adjustment  is  required  throughout 
the  full  humidity  range.  The  capacitive  reactance  is 
associated  with  the  linear  segment  at  the  lowest  humid¬ 
ities  of  the  calibrations  curves  of  Figure  1.  The 
importance  of  this  segment  with  its  relatively  flat  slope 
and  its  shift  to  higher  relative  humidity  ranges  with 
time  becomes  a  problem  as  discussed  later, 

Eddy  correlation  sy stent.  The  fluctuating  wind  vector 
needed  in  (2)  is  measured  with  the  pressure  sphere  anemometer. 
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Fig.  2.  Block  diagram  of  circuit  used  with  barium 
fluoride  film  humidity  sensor  and  its 
associated  thermistor. 
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details  of  which  are  reported  by  Thurtell,  et  al . ,  (1969). 

The  barium  fluoride  humidity  sensor  is  mounted  to  the 
side  of  pressure  sphere  anemometer  and  a  fine-wire 
resistance  thermometer  is  mounted  on  the  other  side 
(Figure  3) .  The  thermometer  provides  air  temperatures 
which,  when  used  in  the  heat  equation  analogous  to  that 
of  (2) ,  gives  the  sensible  heat  flux  density  (Wesely, 
et  al.#  1969).  Figure  4  gives  a  plan  view  of  locations 
of  the  various  components.  In  addition,  c..' d  not  shown 
in  Figures  3  and  4,  a  small  sunshade  was  elevated  15  to 
20  cm  above  the  humidity  sensor.  The  shade  was  used 
since  radiational  heating  often  caused  sensor  temperatures 
to  rise  as  much  as  10  C  above  air  temperature,  which,  ir. 
turn,  caused  the  effective  relative  humidity  of  the  sen¬ 
sor  to  fall  into  its  least  sensitive,  very  dry,  range 
(above  or  near  knee  in  Figure  1) . 

We  anticipated  that  the  spatial  grouping  of  sensors 
was  small  enough  and  the  sensors  had  sufficiently  high 
frequency  responses  to  measure  transport  occurring  in 
small,  high-frequency  eddies  found  within  one  meter  of 
the  surface. 

Calculations  were  made  on-line  by  transmitting  analog 
voltages  of  relative  humidity,  sensor  surface  temperature, 
and  anemometer  pressures  to  an  analog  to  digital  converter 
and  an  Electro-Mechanical  Research  6130  computer  (8001 
Bloomington  Freeway,  Minneapolis,  Minn.  55420)  housed  in 
an  instrument  trailer.  Sampling  rate  was  150  times  per 
second  with  data  acquisition  for  28.5  minutes  of  each 
half  hour  and  summary  data  calculations  and  typewriter 
output  for  the  remainder  of  the  time. 

A  servo-mechanism  rotated  the  instrument  system 
assembly  into  the  wind  as  wind  direction  changed.  Therefore, 


Fig.  3.  (Upper)  Entire  eddy  correlation  system. 


(Lower)  Close  view  of  sensors. 


s*  SA  -  -  • 
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data  was  acquired  automatically,  except  for  sunshade 
adjustment,  gain  adjustments,  and  equipment  maintenance 
which  were  done  during  data  printout. 

Field  trial  site.  During  September  1968,  vapor 
flux  measurements  were  made  at  Hancock,  Wisconsin  over 
snap  beans  (Phaseolus  vulgaris) ,  which  were  30  cm  high. 

A  pressure  cohere  anemometer  with  humidity  sensor  was 
located  60  m  south  of  the  instrument  trailer  with  a  bean 
fetch  of  60  m  to  the  north,  50  m  to  the  east  and  west, 
and  100  m  to  the  south.  A  100  m  to  200  meter  fetch  of 
alfalfa-brome  pasture  extended  beyond  the  beans  to 
shelter  belts  which  were  15  meters  high.  A  second 
anemometer  without  a  humidity  sensor  was  located  10  m 
west  of  the  previously  described  site. 

The  instruments  at  the  humidity  sensor  site  were 
at  an  initial  elevation  of  1.17  m  above  ground  surface. 

This  elevation  was  maintained  until  1030  on  the  14th 
when  it  was  raised  to  2.10  m.  The  2.10  m  elevation  was 
lowered  to  0.75  m  following  0630  on  September  15.  On 
September  20  the  elevation  was  1.17  m.  The  changes  in 
elevations  were  used  to  try  to  detect  the  affect  of  eddy 
size  and  frequency  on  sensor  response. 

Additional  site  instrumentation  provided  two  other 
measurements  of  latent  heat  flux  density  for  comparison 
with  the  eddy  correlation  data.  One  was  evapotrr.nspiration 
measured  with  a  hydraulic  load-cell  lysimeter  [Black,  et  al., 
1969] .  The  other  measurement  was  made  using  the  energy 
balance  equation 


E  =  Rn  -  G  -  Ha  (3) 

The  net  radiation,  R  , 


was  measured  with  a  Funk  net 
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radiometer,  the  soil  heat  flux  density,  G,  was  measured 
using  soil  heat  flux  plates  in  conjunction  with  ther¬ 
mometers,  while  the  sensible  heat  flux  density,  R  ,  was 
obtained  from  the  average  of  the  two  eddy  correlation 
measurements  [Wesely,  et  al,,  1969]. 

Results 

Half-hourly  values  of  latent  heat  flux  density  from 
(2) ,  (3) ,  and  the  lysimeter  are  compared  in  Figure  5  for 
September  12,  13,  14,  and  20.  On  September  20  energy 
balance  data  were  unavailable.  This  figure  also  shows 
half-hourly  mean  values  of  relative  humidity  as  measured 
by  the  barium  fluoride  sensor  and  wind  speed  measured 
by  a  cup  anemometer  mounted  at  1.32  m  approximately  25  m 
southwest  of  the  eddy  correlation  sites. 

Several  conclusions  can  be  drawn  from  these  compar¬ 
ative  data.  First,  there  is  excellent  general  agreement 
both  in  trends  and  in  magnitude  between  the  lysimeter  and 
the  energy  balance  data*  this  confirms  the  validity  of 
our  independent  measurements  used  for  comparisons.  Sec¬ 
ondly,  for  the  most  part,  the  eddy  correlation  data  prior 
to  1000  and  past  1600  on  each  day  show  good  agreement 
with  the  other  two  sets  of  evaporation  data,  while  during 
the  mid-day  period  they  are  one-half  to  one-third  the 
other  data.  Thirdly,  during  data  collection  at  elevations 
of  1.17  in  to  2.1  m,  no  apparent  systematic  differences 
could  be  detected  in  sensor  response  by  comparison  with 
the  independent  methods.  Figure  5  presents  no  data  for 
the  0.75  m  elevation  since  only  two  hours  of  morning 
data  were  collected,  but  these  limited  data  are  also  in 
good  agreement  with  the  independent  measurements. 

Fourthly,  although  wind  speed  is  correlated  with  eddy 


Diurnal  trends  of  latent  heat  flux  density  from 
eddy  correlation  measurements,  from  energy  balance 
and  lysimeter  data,  and  diurnal  trends  of  windspeed 
and  relative  humidity. 
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frequency,  there  was  no  association  between  periods  of 
either  good  or  poor  agreement  and  windspeed.  Finally, 
there  is  a  strong  correlation  between  periods  of  poor 
agreement  and  mean  relative  humidities,  as  seen  by  the 
sensor,  of  less  than  about  45  percent. 

In  Figure  1  the  August  31,  1968  calibration  curve 
shows  a  sharp  break  in  slope  at  42  percent  relative 
humidity.  This  critical  value  was  determined  by  extrap¬ 
olation,  while  the  actual  change  was  most  probably  a 
gradual  one  over  a  ARll  range  of  8  to  10  percent.  The 
slope  of  the  calibration  curve  is  one  of  the  constants 
required  for  the  on-line  computer  program,  and  during 
periods  when  the  mean  relative  humidity  was  near  the 
Xnee  of  the  calibration  curve,  vapor  flux  density  could 
not  be  satisfactorily  computed.  Figure  5  shows  that  all 
periods  of  poor  agreement  occur  when  the  mean  relative 
humidity  was  less  than  or  only  slightly  in  excess  of  the 
critical  value,  42  percent.  A  portion  of  this  decrease 
in  sensor-perceived  relative  humidity  during  mid-day 
periods  is  attributed  to  radiational  heating.  As  wind 
direction  and  sun  angle  varied,  the  small  sunshade  in  a 
fixed  position  relative  to  the  humidity  sensor  frequently 
did  not  shade  the  sensor.  Visual  inspection  of  shade 
orientation  and  adjustment  of  its  position  were  possible 
only  in  the  1.5-min  intervals  at  the  end  of  each  half-hour, 
and  not  during  data  collection. 

The  pressure  sphere  anemometer  measurements  associated 
with  the  humidity  sensor  gave  a  lower  mean  horizontal  wind, 
and  more  negative  o  w'v'  than  those  from  the  anemometer 
without  a  humidity  sensor  attached.  We  doubt  that  this 
was  due  to  spatial  heterogeneity  of  the  row  crop;  it 
most  likely  was  due  to  locating  the  barium  fluoride  element 
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too  far  forward  (see  Figure  4)  where  it  interf erred  with 
the  wind  flow  when  wind  was  from  the  side*  Any  errors 
in  the  cross-wind  measurement  affect  the  wind  coordinate 
transform. 


Recommendations 

The  field  measurements  indicate  that  the  barium 
fluoride  film  humidity  sensor  has  sufficiently  rapid 
response  to  allow  reliable  eddy  correlation  measurements 
of  vapor  flux  within  a  meter  or  less  of  the  surface. 
Modifications  to  the  present  system  should  be:  (1)  The 
sensor  configuration  should  be  changed  from  a  plate  to 
a  cylinder,  with  cooling  tubes  inside  the  cylinder  to 
maintain  the  sensor  at  or  below  ambient  air  temperature. 
Such  a  change  would  permit  temperature  control  of  the 
sensor  so  that  the  operating  point  on  the  calibration 
curve  could  be  kept  away  from  any  "knee" .  The  cylindri¬ 
cal  configuration  also  should  affect  air  flow  around 
the  sphere  less  than  the  plate.  (2)  Since  the  sensor 
calibration  is  altered  by  contamination  with  time,  BaF2 
films  should  be  applied  as  close  as  possible  to  time  of 
use.  The  BaF2  should  be  coated  on  elements  to  which 
thermistor  and  leads  have  been  attached  previously. 

Films  should  be  recalibrated  frequently  during  field 
use.  (3)  Further  tests  should  be  made  of  the  optimum 
sensor  location  with  respect  to  the  pressure  sphere  to 
assure  minimum  interference.  Certainly  the  forward  end 
of  the  sensor  should  be  behind  the  sphere.  (4)  Frequency 
characteristics  of  the  sensor  should  be  established  to 
allow  matching  amplitudes  and  phase  shifts  to  the  wind 
measuring  system. 
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Eddy-correlation  measurements  of  sensible  heat  flux 
close  to  the  earth's  surface  require  fast-responding, 
small  thermometers;  these  can  be  made  with  fine  resist¬ 
ance  wire.  Resistance  thermometers  with  13  u  diameter 
(e.g.  Mcllroy,  1955;  Dyer  and  Maher,  1965;  Hyson,  1968) 
have  better  frequency  response  and  lower  radiation  errors 
tnan  the  commercial ly-available  thermocouples  which  are 
25  u  wire  and  may  have  junctions  larger  than  25  p. .  If 
the  resistance  wire  is  less  than  13  u  in  diameter,  the 
frequency  response  and  radiation  error  is  still  less 
dependent  upon  wind  speed  than  with  13  pt  wire,  and  radi¬ 
ation  error  decreases.  Because  very  fine  wire  thermom¬ 
eters  have  small  radiation  error,  they  can  be  used  to 
measure  average  vertical  temperature  differences  without 
radiation  shielding  or  aspiration;  however,  in  this 
instance  the  rapidly  fluctuating  thermometer  output  is 
filtered  electrically. 

The  purpose  of  this  note  is  to  describe  a  fast,  very 
fine  wire  thermometer  which  is  constructed  easily  and 
which  we  have  found  useful  for  eddy  correlation  measure¬ 
ment  within  0.5  to  1.0  m  of  the  earth's  surface  (Wesely, 
et  al.,  1969) . 

1.  Description  and  construction  of  the  thermometer 

The  supporting  structure  of  the  thermometer,  as 
pictured  in  Figs.  1  and  2,  consists  of  the  frame,  the 
insulating  plug,  and  the  stainless  steel  supporting  tube. 
The  thermometer  element  consists  of  about  55  cm  of 
platinum-coated,  5.6  ^ -diameter,  tungsten  wire^/  wound 


— ^  Sigmund  Cohn,  Mount  Vernon,  N.  Y.  (0.00022-inch 

diameter,  with  about  4  to  7%  weight  platinum  coating) . 


Fig.  1.  Front  view  of  the  resistance 
thermometer. 


Fig.  2.  Resistance  thermometer  details 
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on  the  frame  as  shown. 

a .  Thermometer  frame 

The  frame  consists  of  three  0.16  mm  (0.0063  inch) 

2/ 

enameled  Karma—  wires  which  were  epoxied  to  two  rings 
constructed  from  0.66  mm  (0.0253  inch)  stainless  steel 
wire.  These  rings  are  1.9  cm  in  diameter  and  spaced 
3.2  cm.  The  resistance  wire  is  soldered  to  two  inward 
extensions  of  the  rings.  The  rings  and  their  extensions 
serve  as  electrical  connections  and  provide  mechanical 
support. 

The  rings  are  formed  on  a  mandrel  and  spot-welded. 

The  extensions  are  spot-welded  to  the  rings  and  ail  the 

spot-weld  joints  are  hard  soldered.  The  extensions  are 
3/ 

pretinned—  ,  which  requires  acid  flux;  when  once  pre¬ 
tinned  future  soldering  operations  may  be  done  with 
neutral  and  rosin  fluxes.  Following  pretinning,  the 
stainless  steel  and  any  parts  exposed  to  acid  flux  f nines 
must  be  thoroughly  cleansed  with  soapy  water  and  rinsed 
in  distilled  water. 


2/ 

—  Driver  Harris  Co.,  Harrison,  N.  J. 

3/ 

-  Eutectic  Welding  Alloys  Corporation,  Flushing,  N.  Y„, 
EutecRod  157;  All-State  Welding  Alloys  Co  ,  White 
Plains,  N.  Y. ,  #430  solder. 


66 


b.  Frame  support 

The  frame  support  consists  of  the  stainless  steel 
tube  (9,5  mm  O.D.  x  1,58  nun  wall)  and  the  insulating 
plug,  which  is  held  in  the  tube  with  a  set  screw.  Two 

3 

0.66  mm  diameter  stainless  steel  leads  are  pretinned11/ , 
washed,  and  epoxied  in  holes  drilled  in  the  plug.  The 
stainless  wires  extend  inside  the  tube,  where  they  are 
soldered  to  copper  leads  and  covered  by  heat- shrink 
tubing . 

c.  Winding  the  resistance  wire 

Tungsten  was  chosen  as  a  thermometer  wire  because 
of  its  high  tensile  strength,  but  due  to  the  small  wire 
diameter,  a  load  of  only  seven  grams  will  break  it.  The 
wire  is  best  seen  against  a  dark  background  with  proper 
lighting.  When  the  resistance  wire  is  wound  on  the 
frame,  the  frame  is  attached  with  clips  to  the  end  of  a 
threaded  arbor  which  has  a  pitch  of  2.5  mm  and  which  is 
fixed  in  a  threaded  nut.  About  65  cm  of  the  resistance 
wire  is  unspooled  and  cut  with  masking  tape  folded  to 
the  ends  so  it  may  be  held.  One  end  is  soldered  to  one 
of  the  pretinned  frame  extensions  and  the  other  end, 
weighted  with  the  masking  tape,  hangs  free.  As  the 
arbor  i3  turned  the  resistance  wire  is  pulled  through 
a  stationary  feed  (needle  with  eye  enclosure  cut  away) 
and  wrapped  around  the  frame,  automatically  spacing 
the  windings  at  least  2.5  mm  apart.  Closer  spacing 
may  cause  adjacent  resistance  wires  to  touch  if  the 
Karma  wire  flexes  slightly.  When  ten  windings  are 
on  the  frame,  the  free  end  of  the  wire  is  soldered  to 
the  second  extension.  Before  the  wire  is  wound  on  the 
frame,  the  Karma  cross -struts  are  coated  either  with 
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epoxy  or  with  a  silicone-base  contact  cement—^  to  pre¬ 
vent  the  tungsten  wire  from  sliding  on  the  strut. 

The  platinum-coated  tungsten  resistance  wire  is 
soldered  with  a  sonic  soldering  iron  without  flux  to  the 
pretinned  frame  extension  using  indium  solder-^/.  If  a 
sonic  soldering  iron  is  not  available,  either  a  cut-acid, 
zinc  chloride  flux  or  an  All-State  neutral  420  flux  can 
be  used  with  the  indium  solder*  however,  the  fluxless 
joint  made  with  a  sonic  iron  is  preferable  sincie  no 

T 

electrolytes  are  introduced.  Satisfactory  solder  joints 
cannot  be  made  with  tungsten  wire  that  is  not  coated 
with  solderable  metals. 

To  prevent  misalignment  of  the  fragile  assembly 
after  winding  the  wire,  the  frame  should  be  attached 
immediately  to  the  frame  support. 

: 

f 

2.  Frequency  response  and  radiation  heating 

Chao  and  Sandborn  (1964)  show  that  a  resistance 
wire  thermometer  responds  to  a  temperature  change  as  a 
first-order  system.  For  first-order  systems  the  ampli¬ 
tude  ratios  and  phase  shift  angles,  Q ,  with  sinusoidally 
fluctuating  air  temperatures  are:  | 

I 

A/A  =  (1  +  uiV)-*5  (1) 

r 

A/Aq  =  cos 9  (2) 


4/ 

—  Mystik  Tape,  Inc.,  1700  Winnetka  Ave.,  Northfield, 
Ill.  (Type  A-117) 

^  Indium  Corporation  of  America,  Utica,  N.  Y.,  Indalloy 
solder  #4,  indium  metal. 
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where  r  is  the  time  constant,  u)  is  angular  frequency,  and 
AQ  is  the  output  amplitude  when  u>  =  0. 

Chao  and  Sandborn  derive  an  expression  for  the  time 
constant  of  a  resistance  wire,  neglecting  radiation 
exchange . 


t"1  =  (k/0C)  (tt/L)2  +  (4/D2)  (hD)/pC) 


-  (X/oC)  (4/nD2)  2  I2 


(3) 


where  k  is  the  thermal  conductivity  of  the  wi.re,  I  the 
current  through  the  wire,  L  and  D  are  the  wire  length 
and  diameter,  o  and  C  are  density  and  specific  heat  of 
the  wire,  X  is  the  wire  resistivity,  and  h  «  (k  /D)N 
where  k&  is  the  thermal  conductivity  of  air.  The 
Nusselt  number,  Nu,  for  air  in  transverse  flow  can  be 
found  from 


Nu  =  0.3  +  0 . 51Re 


h 


(4) 


where  Re  is  the  Reynolds  number .  According  to  Grant  and 
Kronauer  (1962) ,  the  Nusselt  number  for  our  extremely 
fine,  long  wire  would  be  slightly  less  than  0.3  in  still 
air;  however  (4)  is  a  sufficiently  good  approximation 
to  calculate  performance  for  field  experiments. 

The  first  term  in  the  right  side  of  (3)  represents 
internal  conduction  along  the  wire  to  the  supports  and 
is  negligible  for  long  wires.  The  second  term  is  pro¬ 
portional  to  the  convection  from  the  wire  per  unit 
temperature  difference  and  is  much  larger  than  the  last 
term,  which  indicates  how  the  temperature  of  the  wire 


effects  the  time  constant.  When  all  of  the  values  for 
our  resistance  wire  are  inserted  into  (3) ,  it  simplifies 
to: 

t  =  1/(1530  Nu  +  5.8)  (5) 

For  a  given  wind  velocity,  (5)  can  be  used  to  calcu¬ 
late  t  for  the  resistance  wire;  then  (1)  can  be  used  to 
find  the  ratio  of  the  amplitudes  for  a  given  frequency 
and  (2)  will  give  the  phase  shift.  For  the  resistance 
wire  used  in  the  thermometer,  the  time  constant  was  cal¬ 
culated  to  be  about  1.5  msec  in  " still" air  and  0.6  msec 
in  10  m  sec-*1  winds,  in  a  laboratory  relatively  free 
of  air  currents,  the  time  constant  of  the  thermometer  was 
observed  to  be  about  1  msec.  Up  to  a  frequency  of  20  Hz, 
reduction  in  amplitude  should  be  less  than  2%  and  phase 
shift  about  10  deg.  Even  at  50  Hz,  less  than  10%  reduc¬ 
tion  is  amplitude  and  a  phase  shift  less  than  25  deg  is 
expected. 

The  frequency  response  may  be  decreased  and  the 
phase  shift  may  be  increased  by  electrically  filtering 
the  analog  signal  from  the  thermometer  bridge.  This  is 
necessary  to  match  the  response  of  a  wind-measuring 
system  when  eddy  correlation  measurements  are  made  or  to 
average  the  signal  when  temperature  gradients  are  meas¬ 
ured.  An  advantage  of  the  fine  wire  is  that  the  effect 
of  the  wind  speed  upon  r  is  negligible  compared  to  the 
total  phase  shift  and  degraded  frequency  response  needed 
to  match  most  eddy  correlation  wind  systems. 

Solar  heating  of  the  fine  resistance  wire  on  the 
thermometer  must  be  dissipated  by  convective  transport. 
Neglecting  other  sources  of  heating,  the  steady  sr.ate 
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energy  balance  can  be  expressed  as 

R  (1-a  )  DIj  =  hDnL  (T-T  )  (6) 

S  S  a 

where  R  is  the  solar  radiation  and  a  is  the  absorp- 

tivity  for  solar  radiation.  Since  this  equation  is  for 

the  radiation  being  absorbed  uniformly  over  the  entire 

cross-section  of  the  wire  and  no  radiation  losses,  the 

calculated  temperature  difference  will  be  an  overestimate. 

-2 

Using  an  extreme  value  of  1400  W  m  for  Rg  and  0.5 

for  a  temperature  differences  for  the  fine  resistance 

s  -1 
wire  are  about  0.15C  in  still  air,  0.09C  in  0.5  m  sec 

wind,  and  0.05C  in  a  5  m  sec”1  wind.  Thus,  it  is  con¬ 
ceivable  that  radiation  could  cause  the  wire  to  heat  up 
as  much  as  0.1C,  but  a  wind  gust  would  not  change  this 
offset  by  more  than  a  few  hundredths  of  a  degree  Celsius. 
This  change  is  negligible  for  eddy  heat  flux  calculations, 
and  often  is  not  significant  when  mean  air  temperature  dif¬ 
ferences  are  needed,  provided  all  the  thermometers  are 
exposed  equally  to  radiation. 

3.  Bridge  design  criteria 

The  temperature  coefficient  of  resistance  of  the 
wire  was  determined  by  measuring  the  thermometer  resis¬ 
tance  in  a  temperature-controlled  kerosene  bath;  it  was 
found  to  be  0.360%  c"1  at  20C,  and  0.350%  c”1  at  30C. 

The  55  cm  of  resistance  wire  wound  on  each  thermometer 
spool  had  a  total  resistance  of  about  1550  ohms,  which 
increased  with  temperature  at  the  rate  of  about  5.5 
ohms  c”1  at  25C. 

The  thermometer  is  measured  in  a  constant-current 
bridge,  with  the  current  low  enough  for  negligible 


71 


self-heating.  Since  the  Nusselt  number  calculated  by  (5) 
for  5.6  p,  wire  only  doubles  as  the  wind  velocity  changes 
from  still  air  to  a  5  m  sec”1  wind,  the  effect  of  con¬ 
vective  heat  transfer  on  self-heating  is  weak.  In  still 
air,  the  maximum  current  allowable  for  less  than  0.01C 
temperature  rise  is  0.33  ma .  This  value  increases  to 
0.35  ma  in  a  0.5  m  sec"1  wind,  and  to  0.50  ma  in  a  5m 
sec"1  wind.  If  the  current  were  0.30  ma,  the  temperature 
difference  would  be  less  than  0.01C  and  would  change  by 

about  0 . 004C  as  the  wind  changed  from  0.5  m  sec”1  to  5  m 
-1 

sec 

In  our  bridge  we  use  a  16.2  V  mercury  battery  (two 
8.1  V,  TR-236R)  across  a  40  k  ohm  resistance  in  series 
with  the  thermometer.  Fixed  resistors  of  similar  value 
and  a  potentiometer  for  balancing  form  the  other  half  of 
the  bridge.  The  bridge  output  is  0.60  C  mv"1.  When  two 
thermometers  are  used  to  measure  vertical  temperature 
differences  for  Bowen's  ratio  measurements  the  second 
half  of  the  bridge  also  is  a  40  k  ohm  resistance  in 
series  with  the  thermometer;  the  40  k  ohm  resistance  is 
comprised  of  a  fixed  resistor  and  a  potentiometer  to 
obtain  balance.  When  used  for  measuring  vertical  temper¬ 
ature  differences,  the  thermometers  are  mounted  on  a 
stand  which  interchanges  their  position  periodically  to 
obviate  zero  errors  (Sargeant  and  Tanner,  1967) . 
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SENSIBLE  HEAT  FLUX  MEASUREMENTS 
WITH  A  YAW  SPHERE  AND  THERMOMETER 
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1.  Introduction 

A  yaw  sphere,  shown  schematically  in  Pig.  1#  when 
directed  into  the  wind  flow,  generates  a  pressure  between 
the  ports  proportional  to  the  product  of  the  horizontal 
and  vertical  winds.  If  this  pressure  is  measured  with  an 
electrical  pressure  transducer  and  if  the  analog  pressure 
signal  then  is  passed  through  a  high-pass  filter  to  drive 
a  resistance  thermometer  bridge,  the  bridge  output  is  pro¬ 
portional  to  u(w't').  The  sensible  heat  flux  density  can 
be  determined  by  integrating  the  bridge  output  and  divid¬ 
ing  the  mean,  u(w't'),  by  the  mean  wind  speed  measured 
with  a  nearby  cup  anemometer. 

The  objective  of  this  note  is  to  describe  this  analog 
system  for  sensible  heat  flux  measurement  and  to  present 
some  comparisons  with  independent  measurements  of  sensi¬ 
ble  heat  flux  density. 

2 .  Equipment  description 

The  description  of  the  equipment  is  helpful  to  a  dis¬ 
cussion  of  the  theory  and  is  given  first. 

а .  Yaw  sphere,  vane,  and  pressure  system 

The  yaw  sphere  was  made  by  drilling  two  1.59-mm  holes 
off-center  through  a  5-cm  plastic  sphere  so  that  the 
included  angle,  9 ,  between  radius  vectors  to  holes  on  the 
sphere  surface  was  45°.  The  sphere  was  mounted  on  a 

б. 35-mm  O.D.  tubular  stem  inserted  in  the  head  of  a  Gill 
propeller  vane  where  the  propeller  normally  mounts, 

(Fig.  1) .  Two  1.59-mm  I.D.  polyethylene  tubes  were  run 
through  the  stem  and  down  the  center  of  the  hollow, 
rotating  shaft  in  the  Gill  propeller  vane  that  drives 
the  azimuth  potentiometer.  These  polyethylene  tubes  were 
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brought  out  through  the  bottom  of  the  vane  housing  and 
attached  to  a  Datametrics,  Model  511-8  Baroeel  capacitive 
pressure  transducer,  A  1-m  length  of  1.59-mm  I.D.  tubing 
was  required  to  connect  each  sphere  port  to  the  pressure 
transducer.  A  fast  resistance  thermometer  (Wesely,  et  al., 
1969a)  was  mounted  on  vane  head  and  located  at  the  side 
of  sphere  as  described  by  Wesely,  et  al.  (1969b) ,  The 
thermometer  leads  were  run  in  parallel  with  the  pressure 
tubing  out  the  bottom  of  the  vane  housing  to  the  ther¬ 
mometer  bridge. 

Although  Wesely,  et  al.  (1969b)  adjusted  the  fre¬ 
quency  response  of  the  thermometer  bridge  amplifier  to 
match  that  of  their  pres sure- sphere  anemometer,  this 
would  have  complicated  our  simple  analog  system;  accord¬ 
ingly  the  yaw-sphere  and  thermometer  have  different  phase 
and  frequency  response.  The  frequency  response  of  the 
yaw  sphere  is  indicated  in  Table  1.  The  thermometer 
relative  amplitude  at  20  Hz  is  about  0.98  with  1  m  sec^1 
winds  and  the  phase  shift  is  about  10°.  Details  of  the 
pressure  transducer  system  and  frequency  response  cali¬ 
bration  methods  can  be  found  in  Thurteli.  et  al.  (1969) . 

Table  1.  Frequency  response  of  the  yaw  sphere,  tubing 
and  Baroeel. 

Frequency,  Hz;  2  4  10  15  20 

Relative  amplitude:  1.0  1,0  0.83  0.63  G.46 

Phase  shift,  degrees *  5  15  35  50  60 

b .  Electronics  and  recording 

The  output  from  the  Baroeel  and  its  signal  conditioner 
is  the  electrical  analog  of  the  pressure,  AP,  between  the 
yaw  sphere  ports.  This  signal  is  passed  through  a  high- 


pass  filter  with  a  17-min  time  constant  end  unity  gain. 

The  output  of  the  high-pass  filter  is  an  analog  of  (AP-AP) . 
This  signal  drives  the  thermometer  bridge  so  that  the 
output  is  an  analog  of  AT (AP-AP) where  AT  is  the  bridge 
temperature  unbalance.  If  the  bridge  is  set  at  a  null 
temperature  very  different  from  the  air  temperature  so 
that  a  large  unbalanced  offset  appears  in  AT,  the 
pe&k-to-peak  range  of  the  fluctuating  bridge  output  is 
unduly  large  and  may  saturate  the  following  electronics ^ 

Tc  facilitate  balancing  the  bridge,  it  can  be  switched 
from  the  high-pass  filter  to  a  battery.  The  output  of 
the  bridge  is  fed  to  a  differential  amplifier  and  thence 
to  either  an  electronic  integrator  or  a  recorder  with 
ball-and-disc  integrator.  The  integrator  is  proceeded  by 
a  5-sec  low-pass  filter  to  decrease  the  transient  response 
and  dynamic  range  requirements. 

3.  Theory  of  operation 

The  pressure  distribution  at  points  on  a  sphere  in  a 
perfect  fluid  with  irrotational  motion  is  given  by  Lamb 
(1932?  sec.  92)  as 

P  =  Ps  +  (p/2) V2[l  -  (9/4) sin2*]  (1) 

where  P  is  the  static  pressure,  u  is  the  density  of  air, 

V  is  the  air  speed,  and  7  is  the  angle  between  V  and  the 
radius  vector  of  the  point.  Schlichting  (I960,  p.  20) 
shows  that  in  real  fluids,  the  pressure  distribution  is 
that  of  ideal  fluids  for  d  <  S5°.  If  the  yaw  sphere  is 
directed  azimuthly  into  the  wind,  then  the  pressure 
difference  between  the  ports  of  the  yaw  sphere  is 


AP  »  P2  -  Px  =  (90/8)|v|  2(sin2*1-sin2*2) 
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This  holds  for  winds  within  a  vertical  angle  a  =  +  (65°-*j0) 
where  <y  is  the  angle  between  the  wind  vector  and  the 
bisect  of  the  ports  and  e  is  the  included  angle  between 
the  ports.  Assuming  that  the  vane  directs  the  sphere 
into  the  wind,  the  components  of  the  wind  vector  with  respect 
to  the  x, z  plane  formed  by  the  ports  and  the  stem  are 


u  «  |v|  cos<y  (2a) 

v  «  0  (2b) 

w  *=  |  v|  sin*  (2c) 

Since  «  (<v  +  *$&)  and  (f2  *  (or  -  *30) 

A P  =  (9/4) (sin  0)  ouw  (3) 


The  electrical  output  of  the  pressure  transducer  is 

Ep  =  M  (AP)  =  bMpuw 

where  b  =  (9/4  sin9)  and  M  is  the  transducer  constant.  The 
output  of  the  high  pass  filter  is 

Ep  =  Ep  -  =  M(AP  -  AP)  (4a) 

Substituting  (3)  into  (4a) # and  using  Reynold's  notation 

Ep  «*  b  M  o  (uw'  +  wu'  +  u'w'  -  u'w')  (4b) 

Since  AT  *=  St  +  AT'  *  AT  +  T' ,  Where  AT  is  the  mean 
bridge  balance  offset,  the  output  of  the  bridge  is 


79 


Efi  -  BEp  AT  =  BEp(AT  +  T')  (5) 

where  B  is  the  bridge  constant.  The  amplifier  output  is 
then 


E0  -  G  Eb  (6) 

where  G  is  the  amplifier  gain,  When  E(  is  integrated  we 
have  from  (4) ,  (5) ,  and  (6) , 

E^  «=  bGBMo  (u  w't'  +  w  u't'  +  u'w't')  (7) 

Assuming  the  last  two  terms  in  parentheses  in  (7)  are 
negligible  as  compared  with  the  first,  the  product  of  the 
sensible  heat  flux  and  the  mean  wind  is 


uH  m  pcp(u  W'T'  +  w  u't'  +  u#w#T* ) 

-  Cp 


(8) 


If  a  cup  anemometer  is  run  near  the  yaw  sphere-thermometer 
assembly  at  the  same  height  to  find  u. 


H»  c  C(9/4)GBM  8ine]_1(E0/u) 


4,  Measurements 

Sensible  heat  flux  density  measurements  were  made  with 
the  yaw  sphere  and  thermometer  during  three  days  in 
September  1968  at  Hancock,  Wisconsin,  The  yaw  sphere  was 
about  95  cm  above  a  crop  of  snap  beans.  We  integrated  EQ 
with  a  ball-and-disc  on  a  strip  chart  recorder  Which  also 
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gave  a  record  of  EQ.  We  used  (9)  to  find  the  heat  flux 
Where  =  45°  for  our  sphere.  The  mean  wind  speed  was 
measured  with  a  cup  anemometer  mounted  at  about  the  same 
height  and  located  15  m  from  the  yaw  sphere.  We  compared 
the  results  with  eddy  correlation  measurements  of  sensi¬ 
ble  heat  flux  density  made  with  a  three-dimensional 
pressure- sphere  anemometer  in  combination  with  a  fast 
thermometer.  We  also  compared  the  results  with  energy 
balance  measurements.  Where  the  sensible  heat  flux  density 
was  found  by  subtracting  measured  soil  heat  flux  density 
and  evaporation  from  the  net  radiation.  Wesely,  et  al. 
(1969b)  describe  the  site,  the  pressure  sphere  anemometer 
and  thermometer  measurements  and  the  energy  balance 
measurements. 

All  measurements  are  given  in  Fig.  2  for  three 
periods.  The  fluctuation  of  the  energy  balance  measure¬ 
ments  is  mainly  due  to  the  lysimeter,  which  is  not  well 

suited  to  measuring  evaporation  over  periods  as  short  as 

2 

30  min  When  peak  evaporation  is  equivalent  to  250  w/m  . 
Also  any  phase  differences  in  the  three  heat  flux  terms 
can  make  for  large  relative  errors  in  sensible  heat? 
this  is  particularly  evident  around  1800  hours. 

The  yaw  sphere-thermometer  results  generally  are 
higher  than  the  other  two  measurements.  The  difference 
between  yaw  sphere-thermometer  data  and  that  of  the  other 
methods  corresponds  more  nearly  to  a  zero  offset  than 
to  a  proportionality  factor.  A  damaging  zero  offset 
could  easily  arise  since  the  mean  bridge  output  voltage 
was  of  the  order  of  80  to  120  iav  during  periods  of  high 
sensible  heat  flux.  In  the  equipment,  which  we  jerry- 
rigged  hurriedly  for  this  teat,  there  were  two  possible 
sources  of  zero  error:  no  particular  precaution  was 
taken  to  avoid  thermal  emf's  in  the  bridge?  also  in  the 
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constructing  of  the  active,  high-pass  filter,  no  special 
attention  was  given  to  avoiding  small  d-c  components  in 
the  output,  and  any  d-c  across  the  unbalanced  thermometer 
bridge  would  result  in  a  zero  offset.  These  features  of 
the  system  can  be  improved  relatively  easily. 

In  view  of  the  success  of  the  preliminary  tests  of 
this  simple  yaw  sphere- thermometer  system,  we  believe  it 
holds  high  promise  for  routine  measurements  of  sensible 
heat  flux  density  as  close  to  the  ground  as  one-meter 
where  other  eddy  correlation  systems  are  not  suitable. 
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ANEMOCLINOMETER  EQUATIONS 
AND  COMPUTER  PROGRAM 

The  instrument,  shown  in  Figures  1  and  2,  consists 
of  a  spherical  head,  3  cm  in  diameter,  mounted  on  a  support¬ 
ing  shaft.  When  in  use  the  probe  is  fixed  in  the  fluid,  with 
the  shaft  axis  parallel  to  the  direction  of  mean  flow  and 
the  sphere  on  the  upstream  end  of  the  shaft.  In  the  follow¬ 
ing  discussion,  all  coordinate  axes,  planes,  and  the  velocities 
are  referenced  to  the  ports  in  the  anemoclinometer  head  and 
the  anemoclinometer  axis. 

The  static  pressures  developed  between  small  holes  drilled 
in  the  spherical  probe  head  are  measured  (Figure  2) .  The 
pressure  difference  between  the  two  holes  in  the  x-z  plane 
is  measured  and  also  between  the  two  holes  in  the  x-y  plane* 
each  of  these  four  holes  is  drilled  at  a  45°  angle  to  the 
axis  of  the  shaft.  In  addition  to  the  above  pressure  ports 
an  upstream  opening  in  the  spherical  head  leads  into  a  Venturi 
centered  on  the  axis  of  the  shaft.  A  small  pressure  tube  is 
placed  in  this  Venturi,  parallel  to  the  probe  shaft,  with 
its  open  end  in  the  upstream  direction.  The  pressure  differ¬ 
ence  is  measured  between  the  pitot  and  eight  reference  ports 
on  the  surface  of  the  spherical  head,  which  are  located  on 
a  circle  at  an  angle  of  47.5°  to  the  shaft  axis. 

The  instantaneous  velocity  component  along  the  z-axis 
and  normal  to  the  shaft  in  the  x-z  plane  of  the  two  vertical 
ports  (plane  Z,  Figure  2)  is  defined  as  w.  The  velocity 
components  in  the  x-y  plane  of  the  two  horizontal  ports 
(plane  Y,  Figure  2) ,  are  u  and  v,  with  u  as  the  component 
parallel  with  the  shaft  (x-coordinate)  and  v  as  the  component 
normal  to  the  shaft  in  the  x-y  plane  (y  component) .  The 
equations  for  the  velocity  vectors  with  respect  to  the 
anemoclinometer  coordinates  are 


Fig.  Pi.  Spherical  sensing  head  of  anemoclinometer 
showing  pressure  ports. 


Fig.  p2 .  Front  and  cross-section  views  of 

anemoclinometer  head,  y-  and  z-coordinates 
shown  on  front  view. 
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u  =  |  v|  cosF' 

cosG  »  |  v|  cosF'  cosG' 

[la] 

v  «  I  V|  cosF' 

sinG  =  |  V| sinG' 

[2a] 

w  =  |v|  sinF' 

[3a] 

|V|  -  (u2  +  v2 

+  w2)*5 

[4a] 

where  F ,  and  G  are  the  elevation  and  azimuth  angles  pro¬ 
jected  on  the  x,z  and  x,y  planes  and  f'  and  G'  are 
complements  of  the  directional  angles  as  shown  in  Figure  3. 
We  can  find 

F'=  arctan  [tan(F  cosG)]  »  F  cosG  [5a] 

The  approximations  in  [la]  and  [5a]  result  in  less  than 
5%  error  at  angles  equal  to  or  less  than  30°. 
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Fig.  P3.  Definition  of  the  components  u,  v#  and  w 

of  the  total  wind  vector  V-  the  elevation  and  azimuth 
angles,  F '  and  G# ,  Which  are  complements  of  the  direc¬ 
tional  angles;  and  F  and  G,  which  are  the  projections 
of  F'  and  G'  on  the  x,z  and  x,y  planes  respectively. 
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An emocli nometer  Constants  (30- mm  Sphere) 


The  pressure  ?v,  measured  between  the  pitot  and  the 
reference  ports  is  related  to  the  true  dynamic  pressure, 


p0*  as 


Pv/P0  =  a  [6a] 

where  PQ  =  (o/2)V  and  where  the  first  anemoclinometer  con¬ 
stant,  a,  is  close  to  1.015  for  the  3-cm  anemoclinometer. 

When  the  centerline  of  the  uw  ports  (anemo¬ 
clinometer  axis)  is  at  an  angle,  F„  with  respect  to  the 
mean  flow  (G  =  0  in  Figure  3)  a  pressure,  P_,  is  developed 

r 

between  the  ports,  where 

Pp  -  bouw  [7a] 

where  b  is  approximately  constant,  and  is  near  1.70  for 
the  3-cm  anemoclinometer  when  2000  <  Re  <  200,000. 

The  ratio  of  Pp/Pv  changes  linearly  with  angle 
(F  *  20°)  as  shown  in  Figure  4. 

(Fp/Pv)/F  =•  c  =  0.057 /deg  =  3. 266 /rad  r«a] 

where  c  is  the  third  probe  constant. 

Note  that  equations  analogous  to  [7a],  and  [8a], exist 
for  the  uv  ports  upon  rotation  through  an  azimuth  angle, 

G  £  20°,  vrhen  F  =  0 


PG  -  bpuv 


[  Sa] 


(PG/PV)/G  -  c  =  0.057/deg  =  3.266/rad 


[  10a] 
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Fig.  P4.  Variation  of  (Pp/Pv)  as  angle  F  is  changed 
with  angle  G  *  0,  or  of  Pq/Pv  »8  angle  G  is  changed 
with  F  =  Cc 


The  above  F,  G  angle  relations  hold  for  the  linear  ^gion 
of  Figure  4  until  the  velocity  vector  is  20°  off  axis. 
Beyond  20°  #  the  relation  is  approximated  by 

F  *  c^(PF/Pv)  -  c'2  [  Ha] 

which  can  also  be  written  as 


U2a] 
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where  F1  =  (Pp/Pv) /c.  Similarly, 

G2  =  c1G1  -  c2  [13a] 

where  G1  =  (PQ/PV)  /c.  In  [12a]  and  [13a],  and  c 2  have 
values  of  =  1,500  and  «  10  degrees  =  0.1745  rad. 
Equations  [12a,  13a]  are  shown  by  the  dashed  line  in 
Figure  4.  Other  approximations  can  be  used,  but  since 
accuracy  at  F^  or  G.?  gr eater  than  30°  is  poor,  there  is 
little  basis  for  choice. 

Wind  Component  Calculations 

To  find  u,  v,  w  using  [la,  2a,  3a],  we  find  V  from 
[6a]  and  find  f'  and  g'  from  calibration  curves  provided 
by  the  Institute  de  Mecanique  des  Fluides  de  Lille  (IMFL) 
for  their  anemoclinometers  (Figure  5) .  These  curves  give 
the  variation  of  (Pp/Pv)  and  (PG/PV)  for  winds  which  are 
outside  of  the  x,  z  and  the  x, y  planes.  IMF^  does  not 
specify  whether  the  symmetric  angles  in  Figure  5  are  the 
angles  F  and  G  projected  on  the  x,z  and  x,y  planes  or  if 
they  are  F*  and  G' ,  which  are  complements  of  the  directional 
angles  (see  Fig.  3) .  We  have  assumed  that  F'  and  G'  were 
the  appropriate  angles-^. 


If  F  and  G  were  the  correct  angles  then  [5a]  would  be 
used  to  determine  the  angles  F#  for  use  in  [la,  2a,  and 
3a] .  We  believe  the  choice  of  f'  and  G'  is  correct 
because  if  we  had  used  F  and  G,  the  values  of  w  and  v 
at  F  =  G  ~  30°  would  have  been  about  10%  smaller  than 
if  F*  and  G '  had  been  used  and  5%  smaller  at  P  =  C  =  20°. 
Our  calculations  of  shear  stress  which  used  F#  and  C*  , 
nc^er  appear  systematically  large  as  would  be  the  case 
if  F  were  the  correct  angle. 
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Fig.  P5 . 
both  F 


Variation  of  (Pp/Py)  and 
and  G  vary. 


<W  as 


The  experimental  data  for  anemoclinometers  shown 
Figure  5  indicates  that  the  (PF/PV)  at  any  angle  G'  i: 
related  to  that  for  g'  =  0  as 

(PF/PV^G'  =  COSl'  *PF/PV*G'=0 


(Wf'  “  coaF'  (pG/PV)  f'=Q 


Similarly, 
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using  [8a]  and  [10a]  we  have 

F'  a  (Pp/Pv)/c  cosG#  [14a] 

G'  =  (PG/Py)/c  cosF*  [15a] 

We  find  F '  and  G '  by  iteration: 

Step  1:  An  angle  F^  is  found  from  [8a], 

Step  2:  An  is  found  from  [12a]  if  >  20°? 

if  £  20°,  F'2  a  F? 

Step  3:  Using  cesF^#  G^  is  found  from  [15a], 
Step  4:  A  is  found  from  [13a]  if  G^  >  20°; 

if  G^  <  20  ,  then  G^  —  G^. 

Step  5:  Using  cosG^#  F^  is  found  from  [14a], 
Step  6:  An  F^  is  found  from  [12a]  if  F^  >  20° ; 
if  F3  <  20°.  . 

Although  further  iterations  could  be  made,  F^  and  are 
within  1%  of  the  values  obtained  by  a  third  loop. 


Azimuth  angle  measurement: 

During  1968,  the  anemoclinometers  were  mounted  on  masts 
which  were  servo-driven  with  slow  motors  to  maintain  orienta¬ 
tion  into  the  wind  with  a  dead-band  of  about  10°.  The  uv 


signal  from  the  anemoclinometers  was  used  for  sensing  direc¬ 
tion.  The  rotation  of  the  masts  was  measured  with  a  poten¬ 
tiometer.  The  azimuth  angle  used  in  the  wind  calculations 
was,  G^,  defined  as 


<g3-g4p> 


[16a] 


where  was  the  angle  of  mast  rotation  measured  by  the 
servo  potentiometer,  G4p  was  the  mean  G^  for  the  previous 
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half-hour,  and  is  the  azimuth  angle  with  respect  to  the 
anemocl inometer  as  found  by  [15a]  in  the  iteration  procedure. 

PROGRAM 


Program  Constants 


Thermometers : 

constants  for  tungsten  wire  air  temperature  thermometer 
B  /  and  bridge  where  T  =  +  B2V4#  Celsius 


constants  for  thermistor  and  the  linearized  bridge  for 
|  measuring  the  surface  temperature  of  the  B  F0  humidity 

aA*  a  * 

sensor.  T„  =  B.  +  B.V^ . 

H  3  4  6 

Heat; 

-3  —i 

CT  =  ncp  =  2.9  cal  cm  K  ,  heat  capacity  of  air 
Vapor  pressure: 

CQ  =  6.108mb  =s  saturation  vapor  pressure  at  zero  Celsius 


C  =  7.5 


constants  in  Teten's  formula  for  cal- 
j  culating  saturation  vapor  pressure,  S, 

C2  =  237.3  Celsius  corresponding  to  a  Celsius  temperature, 

S  (defined  here  S/C  =  10^C1T/^C2+ 

and  on  P14)  0 

=  antiloglQ[C1T/(C2+  T) ] 

Cv  =  slope  of  the  BaF2  humidity  sensor  calibration  curve 
for  operating  range  [ (^relative  humidity) /t volts) 

-3  -1  -1 

Cp  =  4620mb  cm  K  gm  ,  the  specific  gas  constant  for 
water  vapor. 


Wind  and  stress; 

=  1.015  =  a  in  [6a] 

E2  =  3.266rad-1  =  C  in  [8a,  10a] 
E3  -  1.5  a  Cx  in  [12a,  13a] 

Ea  -  0.1745rad  =  C2  in  [12a,  13a] 
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Ee 

E6 

M 


R 


J 


Potentiometer 


G 


3 


constants  to  give 
=  E5V7  +  E6 


(see  [10a]) 

[17a] 


=  range  constant  of  pressure  transducer  which  converts 
output  to  pressure 
=  1.2gm  cm^  =  density  of  air 


Channels 


Channel 

1 

2 

3 

4 


5 

6 


7 


Signal 

Sensor 

Variable 

vi 

Pressure  transducer  #1 

Py  =  (ao/2)  V 

V2 

"  "  #2 

Pp  **  bouw 

V, 

w 

"  "  #3 

PQ  =  bouv 

V4 

Resistance  wire  ther¬ 
mometer  for  air  temper¬ 
ature 

T 

V5 

Ba  rium  fluoride  humidity 
sensor 

H  =  relative 
humidity 

V6 

Ba  rium  fluoride  sensor 
temperature  from  a  ther¬ 
mistor 

th 

V7 

Probe  angle  from  potenti¬ 
ometer  (see  [  16a ,  17a] 

G3 

Initialization  Program 

The  electronics  of  the  pressure  sensors  have  an  electrical 
zerc^VQ,  and  a  full-scale,  Vp,  voltage  readout  for  any  pres¬ 
sure  range.  The  measured  *>e  normalized  to 

(VF~  Fq) .  In  addition,  a  /tight  chamber  is  placed  over  the 
anemoclinometer  sphere  to  shut  out  the  wind  and  short  all 
the  ports  hydraulically;  any  residual  signal,  Vs  on  any 
pressure  range  are  due  to  sensor  offset  and  this  must  be 
accounted  for.  Thus  we  have  a  normalized  voltage  from  any 
transducer 
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Vn  =  (v~  Vs)/(V  V  Uea] 

During  the  initialization  program,  the  V_,  V0  are 

U  r  o 

read  on  each  channel  for  1  to  3  minutes,  averaged,  and 
stored  as  constants  in  the  machine  so  that  the  normalized 
voltages  may  be  calculated. 


Combined  Constants  Used  in  Program 


=  VQ1^  :  M1  ran9e  constant  of  pressure 

transducer  #1  t«„  convert  Vr1  =  (V1-  vsl)/(Vpl~  VQ1)  to 
py(see  [22a]) 

=  M2//^Vf2“  V02^  :  similar  to  but  for  PF 
=  M3/ (Vp3~  VQ3)  :  Similar  to  D1  but  for  PG 

“  D2^1h2 

Note  that  (02/D1e2)  C(V2-  Vs2)/(V1-  vs1)]  ~  (PF/Py)/E2  =  Fx 
=  D3/D^E2  which,  analogous  to  A^,  is  used  to  find  G  from 

<V  VS2)/(V  ''SI5* 

=  (REl/2)  A»x 

Note  that  (V^-  VS1)/A3  -  =  u2  +  v2  +  w2 


V3-  VS3  "  X3 


P1  "  A1«W 


Uc] 

I*]  y 


If  |FX|  >  0.349rad  (20°),  then 


F2|  =  b3|  Fx|  -  E4 


[3A] 


If  I  FjJ  2S  0.349rad  then 


f2I  -  I  pii 


[3B] 


Sign  of  F2  is  the  seme  as  the  sign  of  F^ 

c}  -  a2(x3A1) 


If  |GjJ  >  0.349rad,  then 


e2l  ■  e3I  Gil  -  e4 


[SA] 


If  |g1|  «  0.349rad,  then 


[SB] 


Sign  of  G2  is  the  same  as  the  sign  of  G^. 

If  G2  exceeds  0.69rad  (40°)  set  G2  =  0.69rad  and  record 
number  of  times  G2  >  0.69rad. 


f3  -  a^Xj/Xj.) 


If  |f,I  >  0.349rad, 


P4|  -  e3|  f3|  -  e4 


[7A] 


— ^  Eqs.  [2]  through  [7]  are  from  equations  on  page  f  3 . 
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If  |F3|  «:  0.349rad 


|F4I  “  |P3!  [7B] 

Sign  of  F4  is  the  same  as  the  sign  of  F3.  If  F4  exceeds 
0.69rad,  set  F4  ®  0.69rad  and  record  number  of  time^ 

F4  >  0.69  rad. 


G3  .  e5v7  +  e6 

[8A] 

G4  *  G2  +  <V  G4P> 

[8B] 

where  G4p 
G3  =  mean 

=  +  G3  for  the  previous  run 

G3)  . 

(G2  =  mean  G2, 

L  L 

sinF4  a  w 

[9A] 

V 

cosF4  cosG4  =  A^  u 

[9B] 

sinG4  »  A^  v 

[9C] 

TH  *  B3  +  B4V6 

[10] 

S/Cn  a  10'ClTH/(TH*f  c2^ 

[11] 

4/ 

— ^  NOTE:  in  the  1907  program,  since  the  mast  was  not 
rotated  through  G^,  we  used  G2  in  place  of  G4  in 
[9A,  9B,  9C]. 

See  page  P20  for  alternate  program  equations 
4/ 

—  NOTE:  Teten's  equation  (continued  4/  page  P14) 
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Accumulate  sums  and  calculate  averages  of : 


From 


N1 

= 

(l/n)T  w 

A15 

A3 

[9A] 

N2 

s: 

( 1/n)  T.  u 

A*5 

A3 

[9B] 

N3 

w 

( 1/n)  t  v 

[9C] 

N4 

ss 

(l/n)^.(w 

4 2 

[9A] 

N5 

s 

(l/n)?:(u 

4 2 

[9B] 

N6 

= 

(1/n) E  (v 

4 2 

[9C] 

N7 

s 

(l/n)E  (w 

A3)  (U 

A3> 

[9A],[9B] 

N8 

= 

( 1/n)  ?.  (w 

A3)  (v  A^) 

[9A],[9C] 

N9 

rs 

(l/n)E(u 

A3)  (v  Aj) 

[9B],[9C] 

N10 

= 

<l/n)E[(u  A?*)2 

+  (v  A^)2]*5 

N5.N6 

— 7  NOTE  (Cont.) .  Teten's  equation  is  also  written 

S/eQ  *  antilog  [C^/^  +  C2)  ] 

S  is  the  saturation  vapor  pressure  corresponding  to  the 
BaF2  humidity  sensor  Celsius  temperature,  .  CQ  is  the 
saturation  vapor  pressure  at  OC,  and  Cj,  and  C2  are  con¬ 
stants,  as  given  on  page  P9.  £Tetens,  O.  1930*  "Jber 
einige  meteorologische  Begriffe.  Z.  Geophys.  6:297-309.] 
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Nil  = 

•  (l/n)^(P4) 

[7A],[7B] 

N12  = 

!  (l/n)t(F4)2 

[7A],[7B] 

N13  = 

(l/n)T(G2) 

[5] 

N14  = 

(1/n)  7  (g2)  2 

[5] 

N15  = 

(l/njr^) 

[1A] 

N16  = 

(l/n)>:(X2) 

[IB] 

N17  = 

(1/n)  7  (X3) 

[1C] 

N20  = 

(1/n)  7  (V4) 

Channel  4 

N21  = 

(l/n)7(V4)2 

Channel  4 

N22  = 

(1/n)  7  (V4)  (w  Aj) 

[ 9A] , Channel  4 

N23  =■ 

(l/n)T(v4)  (u  aJ) 

[9B], Channel  4 

N24  = 

(l/n)7(V4)  (v  A^) 

[9C], Channel  4 

K25  = 

(l/n)v(S)/(c0) 

[11] 

N25  = 

d/n)v[  (S)/(co)]2 

1 — » 

V-* 

I-* 

N27  = 

(l/n)7(S/CQ)  (w  A^) 

(v5) 

[11].[9A], 
Channel  5 

N2B  = 

(1/n) 7  (Q/C0)  (u  Aj) 

(V5> 

[  11]  *[9B]  , 
Channel  5 

N29  = 

(1/n)  7  (S/CQ)  (v  A^) 

(v5) 

[11], [9C], 
Channel  5 

N30  = 

(1/n)  7  (v6) 

Channel  6 

H31  = 

(1/n)  7  (Vg)  2 

Channel  6 

N32  3 

(l/n)7  (  :3) 

[3A] 

=  (l/r.)\ (g3) 
=  (l/n)7(G4) 


(l/r)r(G4)‘e 

(l/n)t:(X-)  (VJ 

(l/n)t  (w  A^)  (u  Aj)  (V4) 

(l/n)r(V5) 

(l/n)^(V5)2 

(l/n)^(S/C0)  (Vg) 

(l/n)r[  (S/CQ)  (V5)]2 


[8A] 

[SB] 

[SB] 

[IB], Channel  4 


[9A],[9B] 
Channel  4 

Channel  5 
Channel  5 


[11], Channel  5 
[11], Channel  5 


Output  Calculations  (averages  and  standard  deviations)  ; 


Para- 

meter 

P4  = 

G2  =  G. 

G3  *  G- 

G4  « 

U  =  u 

V  =  v 

W  =r  W 


(N2)/(A^) 
(N3)/(A^) 
(Nl)  /  (Ag) 


TA  =  T  (Bx)  +  (B2)  (N20) 


Stand,  dev. 

[  (N12)  -  (Nil)2]*5 
[  (N14)  -  (N13)2]*5 
[  (N33)  -  (N32)2]®* 

[  (N35)  -  (N34)  2]*5 
{(N5)/(A3)  -  [  (N2)  /  (Ag)  ]  2}  ^ 
{  (N6)  /  (A3)  -  [  (N3)  /  (A^)  ]  2]  ** 
({N4)/(A3)  -  [  (ND/fAij)]2)15 
{  (B^)  2  +  2(B1)  (B2)  (N20) 

+  (B2)2(N21) 

-  [(^)  +  <B2)  (N20)  ]2]  ^ 
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S 

1‘Ti 


(M2  5}  (C0) 

=  l'H  (S3)  +  (B4)  (N30) 

N38 

=  e,mb  [(K40)- 

(N25)  (NSBjKQy,)^) 


{  (M2 6)  (C0)  2_  C  (5I25)(C0)  l2}** 
t  (B3)2  +  2<B3)  (B4)  (N30) 

+  (B4)2(N31) 

-  [  (B3)  +  (B4)  (N30)]V 
[(N39)  -  (N38) 

(CV)[(N41)  -  (N40)23>5  5/ 


5  / 

NOTE:  The  relative  humidity,  H,  is  linear  with 
the  logarithm  of  the  resistance.  The  electronics 
produces  a  voltage  linear  with  H. 

Bx-  «2  -  H'  -  -  klogP.'  =  cvv^ 
v5  “  *5  +  V5 

5=5+  s'  as  found  from  [ill 

s vZ  =  s  v7  +  PvT 

S  5  5 

s"'Vg  =  S  V5  -  svj 

C^,  S'Vj  =  S'H'  =  Cy  ( SV^  -  S  VjT) 

Since  the  vapor  pressure  E  =  Sh 

E  •  SH  =  S  H  +  S'h' 

e  =  s  h  +  c^t  (s/c0W5  ”  (S7C^)  V^] 

F  =  S  H  +  (C^t  (N40)  -  (N25)  (N38)]} 

To  find  E,  S  is  found  in  the  above  output;  H  is  found 
manually  from  V5  given  in  the  output,  using  the  calibration 
curve  for  the  particular  sensor;  the  term  in  braces  is 
calculated  above  as  on  average. 

<E')2  =  Cy[  (Ml)  -  (N40)2) 


•■wvjwasrfrfiL.x 
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[  (60)  (B2)][ (N36)-(N16)  (N20) j  SEE  NOTE  7/ 
J2  [  (60)  (B2)  (Ct)  /  ( A3)  ]  [  (N37)  -  (N7)  (N20)  ] 


—  NOTE:  The  composition  of  the  EU,  EV,  and  EW  term? 

are  similar.  EW  =  Xq'w;  is  the  vertical  vapor  flux 

expressed  as  latent  heat,  where  (60)  is  included  to 
-2  - 1 

give  cal  cm  min  .  In  this  equation 

X  =  [597  ~  0.57T],  T  in  Celsius 

to  give  cal/gm  (latent  heat) 

q  «  e/RvTK*  where  TR  is  Kelvin 

and  is  the  specific  gas 
constant,  so  that 

q  =  e/[Cj;.(T  +  273.2],  T  in  Celsius 


—  NOTE:  and  J2  calculate  the  heat  flux  term 

(oc^)[u  w't'  +  w  u*T'  +  u'w't']  found  by  the  yaw 

sphere  and  thermometer  (Chapter  3) .  is  based  on 

values  which  are  not  corrected  for  pressure  decreases 
with  off-axi*<  winds,  whereas  is  corrected. 
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COORDINATE  SYSTEMS  AND  ROTATIONS 


In  flow  problems,  the  components  of  the  flow  velocity 
must  be  measured  relative  to  some  well-defined  system  of 
coordinates.  In  the  atmosphere,  unlike  duct  flow  Where 
the  direction  of  mean  flow  is  woll  defined,  the  coordinate 
system  used  can  be  selected  arbitrarily.  The  wind  velocity 
components  are  measured  with  respect  to  the  instrument 
coordinates. 

At  the  beginning  of  an  experiment,  the  instrument  may 
be  roughly  aligned  with  respect  to  the  direction  of  mean 
flow;  however,  the  direction  of  the  mean  flow  is  not  known 
precisely  until  the  end  of  the  experiment.  The  instrument 
else  can  be  referenced  with  respect  to  gravity  and  some 
arbitrary  azimuthal  direction. 


Instrument  Coordinate  System 


Let  the  x1#  plane  be  the  plane  of  the  two  horizontal 
ports  of  the  anerooclinometer,  (plane  Y  in  Figure  P2) ,  the 
xl#  Z1  be  the  Plane  of  the  two  vertical  ports  of  the 

anemoclinometer,  (plane  Z  in  Figure  P2)  and  the 
plane  be  the  plane  perpendicular  to  both  the  x^,  y^  plane 
and  the  z 1  plane.  Since  the  relative  positions  of  the 
ports  are  determined  by  precision  machining, the  y^  plane 
ic  probably  perpendicular  to  the  x^#  z^  plane  to  within  a 
ttirute  of  arc,  and  we  can  assume  the  three  planes  are  orthog¬ 
onal. 

The  anemoclinometer  is  aligned  in  the  field  such  that 
the  x^-ayis  is  approximately  the  same  direction  as  the  mean 
wind  vector,  V;  the  y.-axis  is  approximately  normal  to 
gravity  (although  it  could  be  aligned  roughly  parallel  to 
a  sloping  land  or  canopy  surface) ,  the  z^-axis  is  per¬ 
pendicular  to  both  and  2^  to  form  the  third  axis  in  a 
right-handed  coordinate  system,  if  the  anemoclinometer  is 
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moved  in  the  x^,  y^  plane  during  the  measurement  period  to 
keep  the  angle  between  \f  and  x^-axia  within  the  acceptance 
angle  of  the  anemoclinometer  Venturi,  can  redefine  the 
reference  of  the  x^-axis,  for  example  to  be,  true  south, 
or  we  may  also  use  as  a  reference  the  azimuth  of  the  V 
as  determined  during  the  previous  measurement  period.  Wo 
elected  the  latter  in  the  1968  measurements * 

General  Coordina te  Tra n a forms 
New  coordinate  systems  can  be  defined  by  rotations  of 
an  angle  rj  about  the  z^-axis,  the  angle  <5  about  the 
y^-axis,  or  the  angle  r  about  the  x^-axio.  The  angle  r\  is 
positive  as  the  x,y  plane  is  rotated  counterclockwise  as 
viewed  from  the  z^-axis;  the  angle  9  is  positive  as  the 
x,z  plans  is  rotated  clockwise  as  viewed  from  the  positive 
y^-axis;  and  the  angle  ft  is  positive  as  the  z,y  plane  is 
rotated  counterclockwise  as  viewed  from  the  positive  x^ 
axis. 

The  coordinate  system  rotations  are  orthogonal  trans¬ 
formations  and  can  be  represented  in  matrix  form  (e.g. 
Albert,  A.  A.  1949.  Solid  analytic  geometry.  Phoenix  Books 
Science  Series,  University  of  Chicago  Press) .  For  instance, 
if  the  x^,  y^  plane  is  rotated  about  the  z^  axis  by  the 
angle  r\  to  define  new  axes,  x^#  y2#  and 

■  'v  {£)  tiBi 


where  (Z  ) 


/  cosr|  sinr^ 
-  sirni  cost) 
\  0  0 


[IB] 


Similarly,  rotations  through  the  angles  9  or  9  can  be  done 
with  the  orthogonal  matrices,  respectively. 
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and 


<V 


<V 


CO  SB 

0 

sine\ 

0 

1 

0 

sine 

0 

cose  / 

[2] 


[3] 


If  more  than  one  rotation  is  performed,  the  second 
rotation  angle  is  defined  with  respect  to  the  coordinates 
after  the  first  rotation,  A  third  rotation  angle  would 
be  defined  with  respect  to  the  coordinates  after  the  second 
rotation.  Sequential  rotations  are  indicated  by  placing 
the  matrix  of  the  succeeding  transformation  to  the  left  of 
the  matrix  of  the  transformation  that  preceeded  it. 

Natural  wind  coordinate  system: 

A  natural  wind  coordinate  system  may  be  defined  as 
being  a  right-handed  coordinate  system  in  which  the  x-axis 
is  parallel  to  the  mean  flow  with  x  increasing  in  the 
direction  of  the  flow:  thus  w  =  v  =  0,  where  w  and  v  are 
the  mean  wind  components  along  the  z-axis  and  the  y-axis, 
respectively.  The  transformation  from  the  instrument  to 
the  natural  coordinate  system,  requires  the  rotation  through 
angles  r>  and  e  .  Rotation  around  the  x-axis  by  an  angle  B 
will  be  considered  later  (p.R7) ,  but  for  the  present  we 
shall  assume  that  z-axis  is  normal  to  the  land  surface. 

The  instantaneous  wind  components  can  be  separated 
into  mean  values  and  deviations  from  the  means  due  to 
turbulence  and  can  be  represented  in  Reynold's  notation 
as 


R4 


u  b  5  +  u'  [4] 

v  *  v  +  v'  [5] 

w  =  vr  w'  [6] 

and  w  are  time  averages 

where  u,  v/and  u'  ,  v*’  and  w*  are  the  deviations  along  the 
x- ,  y-,  and  z-axes,  respectively.  Since  the  mean  wind 
components  and  v^,  as  measured  in  the  instrument 
coordinates,  usually  are  not  zero,  we  rotate  through  the 
angle  r\  and  then  through  the  angle  e  #  where 


Let 


n  8 

=  arc  tan 

[7] 

4  * 

-  arc  tan 

[V(55 

-2,  Ui 

+ 

CO] 

(CE) 

=  COST) 

s 

v2'*5 

1‘ 

C9] 

(SE) 

=  sirm 

3  ''x 

+  ^ 

CIO] 

(CT) 

=  cose 

a  (u2  +  vj^/tuj 

*<•  V2  + 

cu] 

(ST) 

=  sine 

3 

-2  -2.  »< 

+  vx  +  %#1, 

Then 


[13J 


Therefore, 

u  *  u^CT)  (CE)  +  vx(CT)  (SE)  +  (ST)  [14] 


v  =  v^CE)  -  U1  (SE) 


Us] 
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w  =  w1(CT)  -  (ST) (CE)  -  v1(ST) (SE)  [16] 

Equations  [14],  [15],  and  [IS]  can  be  written  for  the  time- 
averaged  wind  components  or  for  the  fluctuating  components. 


u  =>  (CT)  (CE)  +  v^CT)  (SE)  +  w^ST)  [17] 

u'=  u' (CT)  (CE)  +  (CT) (SE)  +  w^(ST)  [18] 

v'rn  v'(CE)  -  u^  (SE)  [19] 

w'=  w'(CT)  -  u^ (ST) (CE)  -  v'(ST)(SE)  [20] 

By  the  definition  of  9  and  r| ,  v  =  w  =  o. 


By  performing  the  proper  multiplications  and  averag¬ 
ing,  [18],  [19]  and  [20]  can  be  manipulated  to  yield 
following  relationships: 

(\7?  -  (uj)^ (CT)  2  (CE)  2  +  (vJ)2(CT)2(SE)2  +(wJ,)2(ST)2 

+  2u^v'  (CT)  2  (CE)  (SE)  +  2u^w' (CT)  (ST)  (CE)  [21] 

+  2v^w'  (CT)  (ST)  (SE) 

(v')2  -  (v^)2(CE)  2  +  (u^SE)  2  -  2u[v^(CE)(SE)  [22] 

(w')^(q?(CT)2  +  (u^fST)  2  (CE)  2  +(v')2(ST)2(SE)2 
-  2u|w^(CT)  (ST)  (CE)  -  2w^v[(CT)  (ST)  (SE) 

+  2u^v' (ST)  2(CE)  (SE) 


[23] 
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U  w'  =  u;wx  (CE)  [  (CT)  2  -  (ST)2]  -  2u' v|  (CT)  (ST)  (CE)  (SE) 

+  w^(SE)[(CT)2  -  (ST)2]  -  (uj^tCT)  (ST)  (CE)  2  [24] 

-(v^?(CT)  (ST)  (SE)2  +(wp(CT)  (ST) 
u'v'  =  u'v'(CT)[(CE)2  -  (SE)2]  +  ^^[<ST)  (CE) 

-  (ST)  (SE)  -  (uJ)2(CT)  (CE)  (SE)  [25] 

+  (v')2(CT)  (CE)  (SE) 

vw  =  v^w^(CT)  (CE)  -  u'w^(CT)(SE)  -  u^v' (ST)  [  (CE)  2  -  (SE)  2] 
•4-(^(ST)  (CE)  (SE)  -  (v^)2(ST)  (CE)  (SE) 


Similarly,  a  scalar  such  as  temperature  or  water  vapor 
measured  near  the  anemoclinometer  can  be  represented  in 
Reynold's  notation  as  Q  «  q  +  q'  and  covariances  can  be 


corrected  by  the  transform  to  natural  wind  coordinates  as 
follows : 

Q'u'  =  Q'u' (CT) (CE)  +  Q'v|(CT) (SE)  +  Q'w^(ST) 

[27] 

Q'v'  =  Q'vj(CE)  -  Q'u[(SB) 

[28] 

Q'w'  =  Q'w^(CT)  -  Q'u* (ST) (SE)  -  Q'v'(ST)(SE) 

[29] 
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Natural  coordinate  system  with  an  angular  rotation  about 
the  x-axis. 

At  a  site  with  adequate  fetch,  no  divergence,  and 
steady  state  flow,  measurements  indicate  that  in  addition 
to  v  =  w  =  o,  (u')  >  (w' )  1/.  Let.tau  states  that  u'v'  = 

w'v'  =0  2/.  Although  wide  variations  of  (v')2  at  different 
meteorological,  sites  with  similar  conditions  have  been 
reported,  our  measurements  at  a  one-meter  height  indicate 
that  (u'~)~Z  a  (vV  >  (w')2. 

When  u'v'  i  0,  conditions  are  not  ideal;  during  the 
sampling  period  the  horizontal  wind  velocity  tends  to  -in¬ 
crease  as  it  shifts  a  particular  direction  ^/*  When  measure¬ 
ments  indicate  that  u'v'  is  significantly  different  than 
zero,  local  divergence  caused  by  fetch  or  surface  homogeneity 
may  be  occurring,  since  the  coordinate  transform  for  forcing 
u'v'  to  zero  results  in  finite  v  and  w.  A  shift  in  wind 
direction  may  be  accompanied  by  a  change  in  velocity  due  to 
flow  about  large-scale  surface  features;  this  large  scale 
divergence  also  may  affect  u'v'  significantly  over  our 
30-min  sampling  period. 

Since  v  =  0  and  there  is  no  reason  to  expect  v'  to  be 
correlated  with  w'  4/,  measurements 


1/ 

Lumley,  J.  L .  and  H.  A.  Panofsky.  1964.  The  structure 
of  atmospheric  turbulence.  Interscience  Monogr.  Vol. 
12.  John  Wiley  and  Sons,  New  York,  239  p. 

2/ 

Lettau,  II.  H.,  1968.  Three-dimensional  turbulence  in 
unidirectional  mean  flow,  p.127-156.  Iii  Studies  of 
the  effects  of  boundary  modification  in  problems  of 
small  area  meteorology.  U.  S.  Army  Electronics  Command 
Tech.  Rept.  ECOM66-624-A,  156p. 

3/ 

Sutton,  0.  G.  1953.  Micrometeorology  *  McGraw-Hill 
Book  Company,  Inc.  333p. 

4/ 

Sutton,  0.  G.  1948.  Atmospheric  turbulence.  Methuen 
&  Co.  Ltd.  London.  107  pp. 
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of  finite  w'v'  indicate  that  the  z-axis  is  orientated  such 
that  part  of  u'w'  appears  in  w'v'.  By  a  proper  rotation 
of  the  z#y  plane  through  the  angle  q  in  the  natural  wind 
coordinate  system,  w'v#  can  be  set  to  zero,  with  v  and  w 
remaining  zero.  The  result  is  that  the  x,y  plane  is  made 
parallel  to  the  average  slope  of  the  terrain  somewhere 
upwind  from  the  sampling  point. 

Using  (Xq)  from  [3]  and  letting  u^,  v2  and  w^  be  the 
wind  components  after  the  planar  rotation  to  make  w2v2  = 


Therefore, 

u2  =  u 

v2  =  v  (CB)  +  w(SB) 
w2  =  w  (CB)  -  v  (SB) 

where 

CB  =  cosq 

SB  =  sin* 


[30] 


[31] 

[32] 

[33] 

[34] 

[35] 


The  proper  multiplications  and  averaging  of  the  deviation 
of  the  wind  components  results  in 


t 


S* 

ft 


t 


U2W2  =  u'w' <CB)  “  u'v'(SB) 


u'v' 


w'v' 


*  u'v'  (CB)  +  u'w'  (SB) 


=  v'w'[(CB)2  -  (SB)2]  +  (w*)  (CB)  (SB) 


-  (v')2(CB)  (SB) 


and 


(ui)  =  (u') 


=5  u 


v  =  w  =  0 


For  a  scalar  quantity  Q, 


Q'u'  =  Q'u' 


Q'v'  *  Q'v'(CB)  +  Q'w' (SB) 


Q'wi  »  Q'w'  (CB)  -  Q'v'  (SB) 
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(v2'  =  (V')J(CB)2  +  2v'w' (CB)  (SB)  +  (w')2(SB)2  [36] 
(w2)2  =  (w')2(CB)2  -  2w'v' (CB)  (SB)  +  (v')2(SB)2  [37] 


[38] 

[39] 

[40] 


[41] 

[42] 

[43] 


[44] 

[45] 

[46] 


To  make  WjVj  =  0,  we  must  manipulate  [40]  to  get 
(SB)  =  {-  (CB)  +  [(CB)2  +  8(K)2]V4(K)  [47] 


RIO 


where 

K  -  w'v'/t(v')*  -  (w' 7*3  [48] 

The  positive  sign  must  be  used  in  [47]  because  then  (wi)  ^ 

— j — 2  ^ 

is  minimized  and  (v^)  is  maximized,  as  is  desirable  for 
(v^)2  >  (w^)  Since  =  arc  cos  (CB)  is  small,  (CB)  and 
(SB)  can  be  found  in  a  few  iterations  by  first  assuming 
CB  1,  then  solving  [^7]  for  SB,  then 

CB  =*  (K)  /  (SB)  -  2  (K)  (SB)  [49] 

and  repeating  until  sufficient  convergence  obtains. 
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CORRECTED  DATA  FOR  SITE  3.  MAY  2,  1967,  Pagtt  DI0-DI5 


I !  i  VINE  SITE  Mb.  AH  OSU  V4D  HSO  RUW  RUV  RWV  HORU  F  F»U  6  660 

ii  START  HIND  WIND  ST  DEV  REYNOLDS  STRESSES  HIND  ELEV  ANGLE  AZIN  ANGLE 

I:  . •••CM/SEC.* . MM  ...... DYNES/CA2.  ••••••  CM/ SEC  RAD  RAD  RAD  RAD 
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IDENTIFICATION  OF  HEADINGS 
ON  DATA  LISTING 

TIME:  Starting  time.  Pacific  Standard  time  in  1967 

and  Central  Standard  in  1968.  During  1967 , 
the  runs  ended  at  1  minute  and  20  seconds 
before  the  hour  or  half-hour.  During  1968, 
runs  were  for  30  minutes. 

SITE  (1S67) :  The  site  description  is  given  in  Chapter  1, 
as  are  the  instrument  locations  for  April 
26-27,  and  May  2-5.  On  April  22-25,  all 
anemoclinometers  were  at  1  meter  at  the  north 
end  of  the  field. 

SITE  (1968) :  Site  1  was  located  60  meters  south  of  the  in¬ 
strument  trailer  in  a  field  of  snapbeans.  A 
3-cm  anemoclinometer  was  mounted  at  a  height 
of  117  cm.  The  beans  were  25  to  30  cm  high. 

Site  2  was  10  meters  east  of  site  1.  A 
3-cm  anemoclinometer  was  mounted  at  a  height 
of  117  cm  except  following  1030  on  September 
14,  when  the  anemoclinometer  was  moved  to 
210  cm  until  0630  on  September  15.  It  was  at 
75  cm  after  0700  on  September  15  for  the 
remainder  of  the  day. 

The  bean  fetch  was  60  meters  to  the  north, 
50  meters  to  the  east  and  west,  and  100  meters 
to  the  south.  Bevond  the  beans  to  the  south 
was  alfalfa  extending  for  150  meters  to  a  15- 
meter  high  woods.  To  the  west  was  a  100-meter 
alfalfa  field  extending  to  a  10-meter  high 
shelter  belt.  Fetch  to  the  northwest  beyond 
the  beans  was  200  meters  of  low  crops  to  a 
shelter  belt.  To  the  east  was  300  meters  of 
alfalfa  extending  to  a  woods. 
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Mean  wind;  U. 

USD:  Standard  deviation, 

VSD: 

WSD:  "  "  , 

RUW:  Reynold's  stress, 

RUV : 

RWV :  “  ", 

HORIZ.  WIND:  Equivalent  to  anemometer  wind, 

(u^  +  cm/sec. 

P,  ELEV.  ANGLE:  Mean  angle  of  wind  with  x^, 

plane  of  anemoclinometer, 

in  program,  radians.  _  ^ 

FSD,  ELEV.  ANGLE:  Standard  deviation  of  F,  (F* ) ^  . 

G,  AZIM.  ANGLE:  Mean  angle  of  wind  with  the  x^, 

plane  of  the  anemoclinometer, 

{t?2  in  program)  ,  radians. 

GSD,  AZIM.  ANGLE:  Standard  deviation  of  azimuth 

angle,  (G^) 2  in  program. 

WIND  DIR:  Mean  wind  azimuth  direction,  G2  + 

in  program,  measured  clockwise  from 
North,  radians.  (The  listing  is  in¬ 
correct  and  gives  mean  G^;  the  G  AZIM. 
ANGLE,  G2,  should  be  added  to  give 
the  wind  direction) . 

WIND  SHIFT:  Change  in  azimuth  of  mean  direction 

for  one  half-hour  period  from  the 

previous  half-hour,  G4  in  program. 

ETA:  Azimuth  angle  used  in  coordinate  transform, 

arctan  (v^/  u^) ,  radians. 

THETA:  Elevation  angle  used  in  coordinate  trans- 

-  -2  -2 

form,  arctan  [^/(u^  +  v^)  ,  radians. 


(u*  2)  cm/sec. 

(v'2)5*,  cm/sec. 

2 

(w*  )%  cm/sec. 
ou'w'  ,  dynes/cm  . 

.—..I.  ,  1.  2 

ou'v'  ,  dyne s /cm  . 
ow'v'  ,  dynes/cm  . 
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BETA:  Rotation  angle  about  x-axis  (a nemo- 

clinometer  axis) #  to  force  w'v'  =  0, 

see  transform  program  • 

HU:  oc  u't' ,  cal  cm^roin-1  . 

P  ■■  2  M  1 

HV :  pc  v't'#  cal  cm  min“  . 

P  -  _2 

HW|  pCp  w#t'  vertical  heat  flux#  cal  cm  min”  . 

AIR  TEMP*  MEAN:  Mean  air  temperature#  Celsius  • 

— 7  ** 

AIR  TEMP*  ST*  DEV:  (T#  )  . 

EU:  X  u^q7"  #  cal  cm~2min”J'  . 

EV:  X  vTq'  ,  cal  aiT^min*1  , 

EW:  X  w#q'  (latent  heat  of  evaporation) , 

-2  -1 
cal  cm  min 

LIMITS  EXCEEDED  (times  per  100 #000  scans): 

2 

VSQ:  Times  V  voltages  were  negative  and 
set  equal  to  zero  (Program  equation 
CIA]). 

F  :  Times  elevation  angle#  F#  exceeded 

40°  (0.698  rad)  and  was  set  equal  to 
40°  (program  equation  [7])  . 

G  :  Times  azimuth  angle,  G2#  exceeded  40° 
and  was  set  equal  to  40°  (Program 
equation  [ 5] ) . 

In  1967#  the  position  of  the  anemoclinometer 
was  fixed  and  the  azimuth  angle  G  often  was  very 
large*  When  G  was  greater  than  25°,  the  data 
were  discarded*  When  correlation  coefficients 
between  u' #  v* #  and  w'  exceeded  unity,  the  run 
was  discarded* 

The  data  listing  obviously  includes  more 
digits  than  are  experimentally  significant. 
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Notes  on  1968  Data 

No  effort  has  been  made  to  check  the  data  gathered  at 
Hancock,  Wisconsin  during  1968.  The  only  data  excluded 
were  those  Where  notes  indicated  obvious  instrument  failure 
or  when  winds  were  less  than  50cm/sec. 

There  were  times  when  the  azimuth  servo-drive  failed 
and  had  to  be  replaced.  The  accuracy  of  the  azimuth  angle 
may  be  in  doubt  during  preceeding  periods.  Some  notes 
regarding  questionable  periods  are  given  below. 

Site  Is  The  azimuth  potentiometer  was  not 
referenced  during  September 
11-12 

Site  2:  The  servo  system  definitely  mal¬ 

functioned  from  2000h  September  11 
through  0550h  September  12.  The 
motor  required  replacement  at 
1200h  on  September  12.  Prom 
103 5h  September  14  onward,  the 
anemoclinometer  was  on  a  mast  driven 
at  1/3  the  earlier  speed.  We  believe 
performance  was  satisfactory;  how¬ 
ever,  the  slow  response  may  have 
created  larger  error  than  would  be 
observed  with  a  faster  motor. 

General;  At  night  when  winds  were  intermittent 
and  low,  the  uv  signal  to  the  servo 
system  occasionally  was  too  low  to 
actuate  the  motor.  If  wind  direction 
had  shifted  appreciably  following 
a  calm  period,  at  times  the  servo 
system  turned  the 


wrong  diroction  until  it  struck  a 
limit  stop. 


Data  differences  between  site  1  and  site  2  s 

Lower  horizontal  wind,  higher  |pu#w'|  and  more  negative 
eu'v'  generally  are  observed  at  Site  2.  This  possibly  may 
be  due  to  spatial  heterogeneity  of  the  row  crop?  however,  a 
more  likely  possibility  is  that  the  differences  are  due  to 
the  location  of  the  humidity  sensor  (see  Chapter  5) ,  which 
may  have  changed  the  wind  flow  around  the  sphere.  There  may 


be  other  reasons,  not  yet  considered. 
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V 

small,  three-dimensional  pressure-probe  anemometer  (IMFL  anemo- 
clinometer)  was  used  to  measure  the  three  components  of  the  wind 
vector,  shear  stress,  and  the  ratio  of  the  standard  deviation  of 
the  vertical  wind  to  the  friction  velocity  as  influenced  by  atmos¬ 
pheric  stability.  Horizontal  wind  and  shear  stress  have  been  compared 
with  independent  wind  profile  and  shear  stress  meter  measurements. 

The  anemometer  was  coupled  with  a  fa3t  thermometer  for  eddy  corre¬ 
lation  measurements  of  sensible  heat  flux  and  with  a  fast  hygrometer 
for  measurements  of  latent  heat  flux.  The  eddy  correlation  measure¬ 
ments  of  sensible  and  latent  heat  fluxes  were  compared  with  inde¬ 
pendent  energy  balance,  wind  profile,  and  sonic  anemometer-thermometr 
measurements. 
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